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Abstract 



We perform a classification of integrable systems of mixed scalar and vector 
evolution equations with respect to higher symmetries. We consider polynomial 
systems that are homogeneous under a suitable weighting of variables. This paper 
deals with the KdV weighting, the Burgers (or potential KdV or modified KdV) 
weighting, the Ibragimov-Shabat weighting and two unfamiliar weightings. The 
case of other weightings will be studied in a subsequent paper. Making an ansatz 
for undetermined coefficients and using a computer package for solving bilinear 
algebraic systems, wc give the complete lists of 2"^* order systems with a 3^'^ order 
or a 4**^ order symmetry and 3"^^ order systems with a 5*^ order symmetry. For all 
but a few systems in the lists, we show that the system (or, at least a subsystem of 
it) admits either a Lax representation or a linearizing transformation. A thorough 
comparison with recent work of Foursov and Olver is made. 

PACS numbers: 02.30.Ik, 02.30.Jr, 05.45.Yv 
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1 Introduction 



The symmetry approach has been proven to be the most efficient integrabihty test for 
(l+l)-dimensional nonhnear evolution equations [1-10] (see also a recent review [11]). It is 
useful in classifying both scalar evolution equations and evolutionary systems of equations 
(see, e.g. [10]). A mile stone in this direction is the work of Mikhailov, Shabat and Yamilov 
[4,5,7,8] on the classification of 2^^*^ order systems with two components. Their work dealt 
with a large class of systems that are non-polynomial in general and have a nondegenerate 
leading part with respect to x-derivatives. They obtained a complete list of systems 
possessing higher conservation laws, up to some (almost) invertible transformations [7,8]. 
Systems with both higher conservation laws and higher symmetries are believed to be 
integrable by the inverse Scattering method, for short "S-integrable" in the terminology 
of Calogero [12,13]. The aim of this paper is to extend the classification of Mikhailov et 
al. and to make it easily accessible. To be specific, we pursue the following goals with this 
paper: 

• To provide a "user-friendly" complete list of systems without any freedom of non- 
trivial transformations. By that the user does not have to find transformations to 
locate a given system in our list. Trivial scaling parameters are removed. Naturally, 
this is possible only for a much more restricted class of systems than that considered 
by Mikhailov et al. 

• To include systems without higher conservation laws^, but with higher symmetries, 
in the classification. Systems of this sort are believed to be linearizable by an appro- 
priate Change of variables and, if so, said to be "C-integrable" in the terminology 
of Calogero [12,13]. 

• To allow systems to have a degenerate leading part. This means that the coefficient 
matrix of leading terms may have a zero eigenvalue. 

• To classify systems of higher order (3'"'^ order, ...). 

• To classify systems with more than two components. 

Here we mention earlier studies devoted to these extensions, although we do not know 
any work dealing with all these extensions simultaneously. A rather user-friendly list 
of integrable systems of 2"*^ order with two components was presented in [14] (see also 
a similar list in [15]). Some classifications of "C-integrable" systems including coupled 
Burgers-type equations have been reported in [15-18]^. Classification of integrable cou- 
pled KdV-type equations has been performed in [19-22] using the symmetry approach 
and in [23, 24] using the Painleve PDE test. Coupled potential KdV (coupled pKdV) 
equations and coupled modified KdV (coupled mKdV) equations with higher symmetries 
were hsted in [25] (see also [26]). Classification of coupled KdV equations and coupled 
mKdV equations was studied in connection with Jordan algebras in [27,28], where the 

^Here we mean conservation laws that do not depend on x and t explicitly. 

^The list given in [16] seems to be incomplete, because we cannot identify an integrable system of the 
Burgers type (cf. H4.8|l in this paper) with any system in the list. 
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coefficient matrix of leading terms is restricted to the identity. The Painleve PDE test was 
apphed to coupled higher-order nonlinear Schrodinger equations in [29], where intcgrable 
coupled mKdV equations and coupled derivative nonlinear Schrodinger (coupled DNLS) 
equations were obtained. Classification of non-commutative generalizations of integrable 
systems on an associative algebra was addressed in [30] (see also [31-33] for DNLS-type 
systems), while vector generalizations of integrable systems were discussed in [34]. 

In this paper, we investigate evolutionary systems for one scalar unknown u{x,t) 
and one vector unknown U{x,t) = {Ui,U2, ■ ■ ■ ,Un) using the symmetry approach. In 
particular, we classify 2'^'* order and 3"^'* order systems that are polynomial in U and 
their derivatives. This work was initiated by Vladimir Sokolov and the second author 
(T.W.) in [34]. Here, N is an arbitrary positive integer and the product between two 
vectors is defined by the scalar product: 

N 

{d^U, d:U) ^ 5^(9rf/,)(a:C/,), m, n > 0. 

We do not consider constant vectors Cj or matrices Cjk as in Xlj^i C*j(5™f/j) or, for exam- 
ple, ^^k=i^ik{d"''Uj){d^^Uk)- Moreover, we require that the scalar and vector evolution 
equations are truly coupled, that is, U occurs in = ... and u occurs in ^7^ = . . . . 
Classifications described in this paper are restricted to {Xi, X2)- homogeneous systems of 
weight /i. These are systems that admit the one-parameter group of scaling symmetries 

{x, t, u, Uj) — >• {a~^x, a~H, a^^u, a^^Uj), a 7^ 0. 

We consider only systems with Ai, A2 > and a differential order equal to fi. For systems 
with Ai = A2, this would imply the existence of a linear leading part (dispersion), but not 
in the case of mixed systems with Ai 7^ A2. For example, for ^ = 2 and Ai = 2A2, the two 
terms u^x and {U, U^x) have the same weight and a differential order equal to /x. In either 
case, systems having a degenerate leading part are also included in our classification. 

For the scalar case, it was proven in [35] that a A-homogeneous polynomial evolution 
equation with A > and a dispersion term may possess a polynomial higher symmetry 
only if 

A = 2 (KdV weighting), or 

A = 1 (Burgers/pKdV/mKdV weighting), or 

X—\ (Ibragimov-Shabat weighting [36]). 

It was also proven in [35] that any symmetry-integrable^ equation of 2"^^ (3'''^) order in 
the considered classes does possess a symmetry of 3'''^ (5**^) order, respectively. Similar 
results on (Ai, A2)-homogeneous polynomial systems of weight 2 with two components 
were obtained in [14]. Under the conditions of Ai, A2 > 0, |Ai — A2I ^ N>o, a nondegeneracy 
of the hnear part^ and the order of the nonhnear part less than 2, such a system may 
possess polynomial higher symmetries only if Ai = A2 = 2, 1, | or 

^The symmetry-integrable equations are such equations that possess an infinite set of (commuting) 

higher symmetries. 

''For full details, see [14]. 
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'^l — 3, A2 — 3, 

X — 2 \ _ 1 
^1 — 3) ^2 — 3- 

In these classes, any symmetry-integrable system of 2"^^*^ order does possess a symmetry 
of 3^'^ order or 4*^^ order. Since we study (1 + iV)-component systems of 2"^^ and 3'''^ order 
that may have a degenerate leading part, we cannot entirely rely on these results. They 
neither give all possible pairs of (Ai,A2) for integrable cases nor indicate the order of a 
higher symmetry to exist. Nevertheless, in this paper, we concentrate our attention on 
systems that are homogeneous in (Ai, A2) = (2, 2), (1, 1), (i, ^), (i, |) or (|, |) and ji = 2 
or 3. Indeed, as we will see below, there exist a lot of interesting integrable systems in 
these classes. Classifications for other pairs of (Ai, A2) will be reported in a subsequent 
paper. 

The search for integrable systems in this paper is based on the simplest version of the 
symmetry approach [1,2], i.e. the existence of one higher symmetry. It is considered as 
a necessary, but in general not sufficient condition for integrability^. Both "S-integrable" 
and "C-integrable" systems can be detected by the existence of one higher symmetry. 
To do concrete computations using the computer algebra program Crack [40] we as- 
sume the existence of a 3'''^ order or a 4**^ order symmetry for a 2'^'^ order system and 
a 5*'^ order symmetry for a 3''*^ order system. Although the existence of symmetries of 
a specific order may be too restrictive and not necessary for integrability^, it allows us 
to perform exhaustive searches and obtain complete lists for these cases. An overview 
of the performed computations is given in the next section. For the lists generated by 
computer, we first remove unessential parameters by scaling independent and dependent 
variables. We note that a linear change of dependent variables mixing the scalar u with 
components of the vector U would give systems with more than one scalar unknown and 
is therefore not considered. Next, we prove integrability for nearly all listed systems by 
constructing either a Lax representation or a linearizing transformation^. Some systems 
in the lists can be reduced to triangular systems by a nonlinear transformation of depen- 
dent variables. If that is possible then systems contain a closed STibsystem in a nontrivial 
manner. In that case, we first prove the integrability of the subsystem and then discuss 
how to solve the remaining equations. We can reduce the task of proving the integrability 
through establishing relationships among the hsted systems. We construct a rich set of 
Miura-typc transformations, including Miura maps phis potentiation, that connect dif- 
ferent systems in the lists. Thus we have only to investigate one representative for each 
group of connected systems. 

This paper is organized as follows. In section 2, we explain briefly how the hsts of 
systems in this paper are generated by computer. In section 3, we perform a classification 
of 2°'' order and 3"^ order SA'stenis in the Ai = A2 = 2 (KdV weighting) case. The list 

^Some systems of the Bakirov type are known to possess only a finite number of higher symmetries 
[37-39]. However, all such examples are pathological and less interesting, because they are already in 
their given form triangular linear. In this paper, we encounter systems with a higher symmetry that 
are reducible to a triangular form by a nonlinear and non-ultralocal change of variables. It is an open 
question whether such systems are symmetry-integrable in general. 

^As a result, we may miss some integrable cases. 

''In this paper, we are not going to pursue the symmetry integrability. 



7 



of 2°*^ order systems with a 3'"'^ order or a 4**^ order symmetry is empty, while that of 
grd QTi-(^Qj- systems with a 5**^ order symmetry consists of four members. The hst itself is 
already known [34] but we prove their integrability in section 3. 

Section 4 forms the main part of this paper. We classify 2"^^ order and 3'''^ order systems 
in the Ai = A2 = 1 (Burgers/pKdV/mKdV weighting) case. The list of 2°^^ order systems 
consists of three members that generalize the Burgers equation. All these systems possess 
both a 3''*^ order and a 4**^ order symmetry. We can (triangular) linearize two of them 
through an extension of the Hopf-Cole transformation, while integrability*^ of the other 
system, ()4.5|) . remains unproven. We discuss travelling- wave solutions of (a subsystem 
of) this system to indicate its nontrivial nature. The interested reader is referred to 
section l4.2.3l for the details. The list of 3'''^ order systems consists of 25 members, three of 
which are symmetries of the 2^^^^ order systems from the previous list. Consequently, we 
can (triangular) linearize two of the three 3'''^ order systems, while integrability of the other 
grd Qj-f^Qj- system remains to be seen. Another 3'"'^ order system in the list (system ()4.9p ) 
is very close to the latter system and we do not know how to integrate it either. For the 
other 21 (= 25 — 3 — 1) systems of 3'''^ order, we prove that they are integrable or, at least, 
they contain an integrable closed subsystem. We point out that one of the 21 systems is 
the 3'''^ order symmetry of a nontrivial l'^* order system. Miura-type transformations that 
connect 3'''^ order systems in the lists of Ai = A2 = 1 and Ai = A2 = 2 are presented. 

In section 5, we classify 2"^^ order and 3'''^ order systems in the Ai = A2 = | (Ibragimov- 
Shabat weighting) case. The list of 2'^'^ order systems is empty, while that of 3'''^ order 
systems consists of two members. We can linearize both of them through a generalization 
of the linearizing transformation for the Ibragimov-Shabat equation. In section 6, we 
obtain negative results regarding a classification in the case of Ai = |, A2 = |. In section 
7, we perform a classification of 2°^^ order and 3'''^ order systems in the case of Ai = |, A2 = 
|. We obtain one 2'^'^ order system with a 3'''^ order symmetry, two 2°*^ order systems with a 
^th Qj-(^Qj- symmetry and two 3'''^ order systems with a 5**^ order symmetry. Thereby we have 
one 2°*^ order system without a 3'''^ order symmetry, but with a 4**^ order symmetry. All 
the listed systems can be linearized through an ultralocal change of dependent variables. 
Section 8 is devoted to concluding remarks. 

Finally, we would like to mention that our results in section 4 refine and generalize the 
recent work of Foursov and Olver [18,25,26]. Their work focused on polynomial systems of 
two symmetrically coupled nonlinear evolution equations, i.e. symmetric systems for two 
scalar unknowns. They obtained the complete lists of Ai = A2 = 1 homogeneous systems 
of 2^^^^ and 3'''^ order with two higher symmetries of specific orders. Most of the (1 + A^)- 
component systems listed in section 4 generalize two-component systems of Foursov- 
Olver, up to a linear change of dependent variables. To see this, we remark that, because 
of our assumption on the admissible multiplications, the evolution equation for the scalar 
u is even in the vector U, while the equation for U is odd in U. Therefore, we can 
symmetrize our systems in the special case A^ = 1 through the linear change of variables: 
u = a{q + r), U = b{q — r) where a and b are nonzero constants. 

On the basis of this re-formulation, we compare in section 4 our lists with those of 
Foursov-Olver. A brief summary of the comparison results is as follows: 

^We mean the existence of either a Lax representation or a Unearizing transformation. 
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• Any system in Foursov-Olver's lists^ corresponds to the N = 1 case of one or two 
systems in our lists. This means that, after the hnear change of dependent variables 
mentioned above, their two-component systems always admit (1 + iV)-component 
generalization (s) preserving the integrability. This result is quite interesting, but 
unlikely to hold true in general for other classes of two-component systems. 

• Some systems in our lists do not have any counterpart in Foursov-Olver's lists. They 
are systems that become the trivial equation Ut — under the reduction U — 0. 
Such systems were excluded from consideration in the work of Foursov-Olver by 
their assumption on strong nondegeneracy of the linear part (see, e.g. section 2 
of [18]). In this respect, our lists are richer than Foursov-Olver's lists even in the 
N = 1 case. 

Besides extending Foursov-Olver's lists [18,25,26], we prove the integrability of many 
systems in their lists for the first time. We also correct errors in [18,25,26] and point out 
overlooked references in which some systems in their lists were studied earlier. 



^We mean the lists of systems that are not reducible to a triangular form by a linear change of 
dependent variables. 
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2 Computational aspects 



Before describing the classification results in sections OHZl in this section, we would like 
to make some comments on the computations performed. 

As the first step a homogeneous ansatz for a system consisting of a scalar equation ut = 
. . . and a vector equation Ut = ... is generated together with a system of higher symmetry 
equations Ur = . . . , Ur = . . . . Each term has a different undetermined coefficient. We 
assume that these coefficients do not depend on (the number of components of U) , and 
that N is not fixed at any specific value^°. 

Computationally more expensive is the formulation of the symmetry conditions uit^r] = 
0, U[t.T] = 0. For low values of Ai,A2 and high differential order the right-hand sides of 
the system and the symmetry involve many terms, and in addition each of the terms has 
an increasing number of factors. Higher-order x-derivatives of such terms cause a large 
expression swell, too large to compute the commutators in one step. We therefore perform 
the computation of U[t^r] and U[t^r] in stages. Because right-hand sides of the system and 
the symmetry do not involve dt, dr, substitutions of Ut, Ut, Ur, Ur in the commutators are 
done only once. Consequently, commutators are linear in the coefficients of the system 
and coefficients of the symmetry. To exploit this linearity we partition 



where the expressions Fi,Gi, Hi, Ki, P^, Qi contain only terms with a total degree i of all 
scalar vector products of U and x-derivatives of U (for example, {U, Ux)U having degree 
1). By using the observation that the number of scalar vector products in a term does 
not change when a term is differentiated we can compute each Pi independently through 



and similarly for each Qi. Because they are the only terms that have z-th degree powers 
of scalar vector products, all Pi, Qi must vanish identically. After one single Pi or Qi 
is computed, it can be split and some of the consequences, like the vanishing of some 
coefficients, can be used to simplify Fi,Gi,Hi, Ki before computing the next Pj and Qj. 

For large problems (low A and high differential order) the computation of Qi was still 
too memory intensive^^ so that another partitioning of the computation was implemented. 
In this level of partitioning, first those terms in Fi,Gi, Hi, Ki which can contribute to 
the highest derivative vectorial factor d^U in Qi were considered. Let us call the partial 
commutator that comes out of this computation Cij. From Ci,j all the terms with vectorial 
factor diU are extracted, split, and the remaining terms from Cij (with a vectorial factor 

-"^^The arbitrariness of is crucial for functional independence of the scalar products 
(0 < m < n) [34]. 

^"'^By splitting we mean the extraction and setting to zero of all coefhcients of all products of all powers 
of all scalar and vector functions and their derivatives. 

""^^The computation was too memory demanding to be performed on the computers available to one of 
the authors in 2000. 



(2.1) 




j=0 
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of order < j) are carried over for the computation of the next terms with derivative 
vectorial factor di~^U in Qi}^ 

In this way the computation of single large expressions M[t.r], t^[t,r] is avoided and 
replaced by the computation of many partial commutators resulting in bilinear algebraic 
equations for the undetermined coefficients. 

To each list of conditions is attached a list of inequalities which have to be fulfilled 
by any solution. A first inequality results from the requirement that at least one of both 
equations involves at least one x-derivative of the required order (2 or 3). Similarly the 
symmetry equations have to involve at least one x-derivative of the required order and the 
right-hand sides of both symmetry equations must be non-zero. Two further conditions 
prevent the generation of triangular integrable systems by requiring that U occurs in 
Ut = . ■ . and u occurs in Ut = ... . 

The solution of the overdetermined bilinear algebraic systems was accomplished with 
the computer program Crack written for the solution of overdetermined algebraic but 
also differential systems. One technique that proved to be quite useful in general, espe- 
cially for the solution of larger systems with Ai = A2 = 1 and Ai = A2 = |, is an equation 
shortening method described in [41]. 

The following tables give an overview on the complexity of computations. These 
have been performed on a 1.7GHz Pentium 4 PC running the computer algebra system 
REDUCE 3.7 in a 120 MB session under Linux. Quoted execution times are sensitive to 
settings of computing parameters and should be taken only as rough indicators. 

For Ai = A2 = 1 and orders 3+5 the large number of solutions has the consequence that 
the system of algebraic conditions does not simplify so readily and is more complicated to 
solve. Hence the program has more often to impose case distinctions where an unknown 
is assumed to be at first zero and then non-zero. As a result solutions may be found 
which can be unified into a single solution. This is the case if, for example, one solution 
Si includes the condition 017 = while the other solution 5*2 requires 017 7^ for some 
undetermined coefficient a^, and if setting 017 = in S'2 makes both solutions equivalent, 
in the system and in the symmetry. Sometimes a substitution like 017 = in S'2 may 
cause a division by zero which can be avoided by re-parametrizing S'2. An algorithm and 
its implementation in a computer program analyzing such situations have recently been 
developed and applied. 

On the web page 

http : //lie . math . brocku. ca/twolf /htdocs/sv/over .htmll 
one can inspect the original systems of conditions and the solutions as well as download 
them in machine readable form. In addition to investigating systems of differential orders 
2 + 3 (for system + symmetry), 2 + 4 and 3 + 5, we also investigated orders 1 + 2 and 
1 + 3. The main purpose was to recognize whether a 2'^'^ order or a 3'''^ order system is 
actually the symmetry of a nontrivial 1*^* order system. Details can also be found on the 
above-mentioned web page. The package Crack can be obtained from 
http : //lie . math . brocku . ca/ twolf / crack/ . 

-"^^More details can be obtained at request. 
Although the program Crack originally produced 4 solutions, we could easily recognize that one 
solution is a special case of another. 

^^The computation involved one manual interference. 
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Ai, A2 


2,2 


1,1 


1 1 
2 ' 2 


1 2 

3 ' 3 


2 1 

3 ' 3 


^ of unknowns in the system 


5 


10 


15 


10 


13 


7^ of unknowns in the symmetry 


6 


21 


36 


24 


22 


total ^ of unknowns 


11 


31 


51 


34 


35 


^ of conditions 


13 


66 


149 


102 


114 


total 7^ of terms in all conditions 


34 


341 


1093 


529 


694 


average # of terms in a condition 


2.6 


5.2 


7.3 


5.2 


6.1 


time to formulate alg. conditions 


0.5s 


1.8s 


8s 


3.2s 


6.3s 


time to solve conditions 


0.5s 


29s 


29s 


45s 


22s 


7^ of solutions 





3 








1 



Table 1: Computations in the orders 2 + 3 problem for the 5 weightings 



Ai, A2 


2,2 


1,1 


1 1 
2 ' 2 


1 2 
3 ' 3 


2 1 

3' 3 


7^ of unknowns in the system 


5 


10 


15 


10 


13 


^ of unknowns in the symmetry 


12 


39 


79 


54 


66 


total # of unknowns 


17 


49 


94 


64 


79 


# of conditions 


26 


123 


313 


215 


276 


total # of terms in all conditions 


77 


770 


3096 


1462 


2435 


average of terms in a condition 


3.0 


6.3 


9.9 


6.8 


8.8 


time to formulate alg. conditions 


Is 


5s 


48s 


13s 


48s 


time to solve conditions 


0.4s 


lm58s 


3m44s 


lm23s 


3m40s 


# of solutions 





314 








2 



Table 2: Computations in the orders 2 + 4 problem for the 5 weightings 



Ai, A2 


2,2 


1,1 


1 1 

2 ' 2 


1 2 

3' 3 


2 1 

3' 3 


# of unknowns in the system 


6 


21 


36 


24 


22 


^ of unknowns in the symmetry 


17 


74 


164 


115 


126 


total # of unknowns 


23 


95 


200 


139 


148 


# of conditions 


50 


386 


1154 


798 


955 


total # of terms in all conditions 


218 


5000 


27695 


12694 


17385 


average # of terms in a condition 


4.4 


13 


24 


16 


18 


time to formulate alg. conditions 


5s 


2m52s 


2h7m 


23m45s 


41ml8s 


time to solve conditions 


6.5s 


5h47m 


Iday^^ 


lh20m 


lh7m 


# of solutions 


4 


25 


2 





2 



Table 3: Computations in the orders 3 + 5 problem for the 5 weightings 
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After all solutions have been determined the task of proving intcgrability follows. In 
the process of identifying and classifying some of the constant coefficient systems, the 
Mathematica package "InvariantsSymmetries.m" [42] has been used to compute conser- 
vation laws and higher symmetries. 
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3 The case Ai = A2 = 2 — coupled KdV equations — 



In this section, we classify 2^^^^ order and 3^'^ order systems in the Ai = A2 = 2 (KdV 
weighting) case. In the first part (section EH)), we present a complete list of such systems 
with a specific order symmetry (the list is already known, see [34]). In the second part 
(section we prove the integrability of the listed systems. 

3.1 List of systems with a higher symmetry 

The general ansatz for a Ai = A2 = 2 homogeneous evolutionary system of 2^^^^ order for a 
scalar function u and a vector function U takes the form 

(ut2 = aiUxx + a2u'^ + a3{U,U), 
I Ut^ = aiUxx + a^uU. 

The following constraints guarantee the order to be 2 and the system not to be triangular: 

(oi, 04) 7^ (0, 0), 03 7^ 0, as ^ 0. 

Similarly, the general ansatz for a 3'''^ order system takes the form 

r ut^ = biUxxx + b2UUx + bs{U, Ux), 
I Ut.^ = bJJxxx + b^UxU + bQuUx, 

for which the following constraints guarantee the order to be 3 and the system not to be 
triangular: 

(61,64)^(0,0), 637^0, (65,66)7^(0,0). 
However, we relax these constraints as 

(61, 64) ^ (0, 0), (61, 62, 63) ^ (0, 0, 0), (64, 65, 66) ^ (0, 0, 0), 

when we consider a 3'''^ order symmetry for a 2°^^ order system, as stated in section |21 We 
omit the general ansatz for a 4*^^ order or a 5**^ order system here (in the Ai = A2 = 2 
case) and hereafter (in the other weightings) because of its increased length. However, it 
is available on the above-mentioned internet site. 

Proposition 3.1. No 2"'' order system of the form i|) with a 3'"'^ order symmetry of the 
form 1)5'. or a 4*'' order symmetry exists. 



Theorem 3.2. Any 3'"'^ order system of the form with a S*'' order symmetry has to 

coincide with one of the following four systems up to a scaling of t3,x,u,U {we omit the 
subscript oft^): 

u. = {U,Ux), ^3^3^ 

Ut = Uxxx + UxU + 2uUx, 
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Ut = Uxxx + QuU^-Q{U,Ur,), ^^^^ 

Ut = Uxxx + QuxU + 6uUx, 

Ut = Uxxx + ^UUx + 3{U,Ux), , ^. 

(3.5) 

Ut = UxU + uUx, 

Ut = Uxxx + Quux - l2{U,Ux), ^. 

(3.6) 

Ut = —2Uxxx ~ QuUx- 

All systems (j3.3|) - (|3.6p admit the reduction U = 0. From this viewpoint, (j3.3p is a 
generalization of the trivial equation Ut = 0, while (j3.4j) - (|3.6p are generalizations of the 
KdV equation. 

3.2 Integrability of systems ({^-(jHTHH 
3.2.1 System 

System ()3.3|) is a multi-component generalization of one of the Drinfel'd-Sokolov systems 
[43,44]. The integrability of this system has been established in the literature [45-47]. 



3.2.2 System 

System ()3.4|) is known as a Jordan KdV system [27,28,34,48]. Let us briefly summarize 
its integrability. It is well-known that the matrix KdV equation, 

Qt = Qxxx + 'i{Q'^)x, (3-7) 

admits a Lax representation [49-52]. Then, system (j3.4|) is also integrable, because it is 
obtained from (|3.7|) through the following reduction: 

N 

Q = ul + ^ UjCj. 

i=i 

Here 1 is the identity matrix and {ei, . . . , e^r} are mutually anti-commuting matrices that 
satisfy the condition, 

{ei, ej}+ = e^ej + eje^ = -25ijl. 



3.2.3 System (jSISD 

System ()3.5|) is a multi-component generalization of the Zakharov-Ito system [52,53] and 
corresponds to a special case of the coupled KdV equations considered by Kupershmidt 
[54]. Introducing a new variable w by w = A/(^^rf^) we find that ()3.5|1 contains the 
original Zakharov-Ito system, 

( Ut = Uxxx + ^UUx + ^WWx, , ^ 

[ Wt = [uw)x. 
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Therefore, system (13 .Sp is a triangular system that consists of the Zakharov-Ito system 
and the hnear equation for U with Zakharov-Ito-system-dependent coefficients. 

To demonstrate the integrabihty of the whole system ()3.5p . we first summarize a 
Lax representation for the Zakharov-Ito system [52,55,56]. We consider a pair of linear 
equations for a scalar function ip, 

(i^.. = {C + q + C'r)ij, 

= (4C - 2g)^, + q,^, 

where ( is the spectral parameter. Then the compatibility condition ip^^t = iptxx for ()3.9|) 
implies the following system: 

{qt = q-xxx - Qqqx + 'ir^, 

\rt = -Aq^r - 2qr^. 

This system coincides with the Zakharov-Ito system ()3.8|) through the change of depen- 
dent variables, q = —u/2, r = — 3w^/16. It is noteworthy that the quantity l/ip"^ in the 
limit C ^ obeys the same evolution equation as that for U, namely, Ut = {uU)^- 

Next, we fix a solution of subsystem ()3.8|) and discuss solutions of the linear equation 
for U. For the sake of simplicity, we assume that w{x,t) in the fixed solution is not a 
trivial function. Then, noting the relation Wt = {uw)x, we obtain the following solution 
to the equation for U: 

Uj = w-f,{^j wdx'), j = l,2,...,iV. 

Here fi{z), . . . , f^iz) are arbitrary functions of z, except that they must satisfy one 
constraint, Ifji^)]^ = 1) due to the relation {U,U) = w^. For the case in which 

w{x,t) is identically zero, we mention some references in section 14.2.91 



3.2.4 System (HHD 

System ()3.6p is a multi-component generalization^^ of the two-component KdV system 
[ (j3.6|) with = 1] proposed by Hirota and Satsuma [57]. Actually, the Hirota-Satsuma 
system is also understood as an example of the Kac-Moody KdV systems studied indepen- 
dently by Drinfel'd and Sokolov [43,44]. A Lax representation for the Hirota-Satsuma 
system was constructed in [59]^''. Recently, it was generalized to the three-component 
case [ ()3.6|) with A^ = 2] by Wu et al. [61]. Let us demonstrate that ()3.6|) admits a Lax 
representation in the general case of A^-component vector U . We consider a set of linear 
equations for two column-vector functions '0 and 0, 

I 0, = + 2P0,. - 4i?V. - + 2i?,.V. 

^^This multi-component generalization was proposed in [58], but the integrability was not discussed in 
that paper. 

^^Vladimir Sokolov commented that the Lax representation was reported earlier in the Russian paper 
[60], which is not accessible to the authors. 



16 



Here, ( is the spectral parameter and P, Q and R are square matrices with the same 
dimension. The compatibihty conditions il^^xt = '^txx^ 4>xxt = imply a system of 
three matrix equations, 

Pt = Pxxx + KP% - HQR)x, 

Qt = -2Qxxx - QQxP + s[Px, Q], (3.10) 

Rt = —2Rx^^ — 6RxP + 3[Px,R], 
together with three constraints, 

[P,Q] = 0, [P,R] = 0, [Q,R]x = 0. 
If we consider the reduction, 

N-l N-1 

P = ul, Q = f/il + ^ Uj+iBj, R = Uil-'^ Uj+iej, {ei, ej}+ = -2Sijl, 

i=i j=i 

the three constraints are automatically satisfied and system (j3.1Up is reduced to the multi- 
component Hirota-Satsuma system (|3.6p . 
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4 The case Ai = A2 = 1 — coupled Burgers, pKdV, 
mKdV equations — 

In this section, we classify 2°^^ order and 3'''^ order systems in the Ai = A2 = 1 (Burg- 
ers/pKdV/mKdV weighting) case. In the first part f section l4.ip . we present complete 
lists of such systems with a specific order symmetry. In the second part f section |4.2|) . 
we discuss the integrability of the listed systems and compare them with Foursov-Olver's 
two-component systems [18,25,26] through symmetrization as stated in the introduction. 

4.1 Lists of systems with a higher symmetry 

The general ansatz for a Ai = A2 = 1 homogeneous evolutionary system of 2°*^ order for a 
scalar function u and a vector function U takes the form 



The following constraints guarantee the order to be 2 and the system not to be triangular: 



' Mt3 = hiUxxx + h2UUxx + hul + hiU^Ux + hu^ + bQUx{U, U) 
+ h-ju{U, Ux) + 68(f/, Uxx) + h{Ux, Ux) + 6iow'(t/, U) 

Ut^ = buUxxx + bisUxxU + buUxUx + h^uUxx + biQUUxU 
+ biju^Ux + bis{U, U)Ux + bi,{U, Ux)U + b2ou^U 

[ +b2Mu,u)u, 



for which the following constraints guarantee the order to be 3 and the system not to be 
triangular: (61, 612) 7^ (0, 0) and at least one of 63, . . . , 611 and one of 613, ... , bn, 620 ? ^21 
must not vanish. However, when we consider a 3'''^ order symmetry for a 2^^*^ order system, 
we relax these constraints as follows (cf. sectional): (&i,&i2) 7^ (0,0) and at least one of 
61, . . . , 611 and one of 612, ... , &21 must not vanish. 

Theorem 4.1. Any 2"*^ order system of the form ( |^.i| ) with a 3^'^ order symmetry of the 
form l\4-^ has to coincide with one of the following three systems up to a scaling of 
t2,x,u,U {we omit the subscript oft2): 




(4.1) 



(ai, oq) ^ (0, 0), (04, 05) 7^ (0, 0), (ay, ag, ag) 7^ (0, 0, 0). 
Similarly, the general ansatz for a 3''*^ order system takes the form 




(4.3) 




(4.4) 




(4.5) 
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Proposition 4.2. Any 2^^ order system of the form with a 4*'* order symmetry has 

to coincide with one of the three systems i\4-^ - ^T~^ W to a scaling oft2,x,u,U. 

All systems ()4.3|) - ()4.5|) admit the reduction U = 0. From this viewpoint, ()4.3p - ()4.5p are 
considered as generalizations of the Burgers equation. 

Theorem 4.3. Any 3^'^ order system of the form i \4-^ with a 5*'^ order symmetry has to 
coincide with one of the following 25 systems up to a scaling of t3,x,u,U {we omit the 
subscript oft^): 

' ut = a{u^^^ + 3uu^^ + 3ul + 3u^u^) + u^{U, U) + 2u{U, U^) 
+ 2{U,U,,)+2{U,,U,), 

Ut = U^^^ + 1(1 - a)u^^U + lu^U^ + luU^^ + 1(1 - 2a)uu^U (4-6) 
+ \v?U^ - {U, U^)U + |(1 - Aa)u^U - \u{U, U)U, 

a : arbitrary, 

Ut = u^xx + 3uu^^ + 3ul + 3u'^u^ + Ux{U, U) + 2u{U, U^) 

+ 2{U,U,,)+2{U,,U,), (4.7) 

Ut = -\u^xU - \uu^U - (f/, U^)U - \u^U - \u{U, U)U, 

Ut = u^^^ + 3uu^^ + 3ul + 3u^Ux + M^(f/, U) + 2u{U, U^) 

+ {U,U^cc) + {U^,U^), (4.8) 
Ut = \u^xU + u^U^ + uu^U + u'^U^ + i(f/, U)U^ + i(f/, U^)U, 

Ut = u^^^ + 3uu^^ + 3ul + 3u^u,j, + u^{U, U) + 2u{U, U^) 

+ {U,U,,) + {U,,U,), (4.9) 

Ut = \u^xU + M^^f/^ + uu^f/ + u^U^ + (f/, U)U^, 



Ut = 3u,{U, U) + 3(t/, U^^) - 3{U, U)\ 
Ut = U^xx + u^JJ + uJJ^ - 3{U, U^)U, 

Ut = 2uAU, U) + 2(f/, U^,) - (f/,, U^) - 2{U, U)\ 
Ut = U^xx + u,,U + 2u,U, - 2{U, U)U, - 2{U, U,)U, 

Ut = u,{U, U) + 2u{U, U,) + {U, U^,) + (f/,, U^), 

Ut = U^xx + u^^U + uJJ^ - 2uu^U - u^U^ + {U, U)U^ - (f/, U^)U, 

Ut Uxxx ~l~ 2Ux ~l~ '2^\Ux^ Ux) , 
Ut = UxUx, 

Ut = Uxxx + ^ul + 2aUx{U, U) + a(f/, Uxx) + cl{Ux, U^) + h{U, U)^, 
Ut = UxxU + 2uxUx + a(t/, U)Ux + a{U, Ux)U, (a, b) ^ (0, 0), 



(4.10) 



(4.11) 



(4.12) 



(4.13) 



(4.14) 
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(Ut = Uxxx + - 3([/a;, Ux), 

\Ut^ Uxxx + QUxUx, 

j Ut = Uxxx + + Ux{U, U) + Uxx), 

\Ut = Uxxx + 3Ma;a;t/ + ?>UxUx + (t/, 

{ut = Uxxx + 3?x^ + 2ux{U, U) + ?7,,) + \{Ux, Ux), 
\Ut = Uxxx + QuxxU + 6M,,[/a; + 2([/, 

r = ^ixx. + ^ul + 4ii,(C/, [/) + 2(C/, + + §([/, C/)2, 

I C/t = -2C/^^^ - %UxxU - QuxUx - 4(C/, Ux)U, 

{ut^ Uxxx + ul- 12(C/, Uxx) + 12([/., - 4([/, 

I C/t = 4C/,,, + UxxU + 2ii,C/, + 4(C/, + 4(C/, [/,)[/, 

\ut = Uxxx - lu^u^ + lu^{U, U) + C/^) + ([/, [4a;) + ([/^, Ux), 

\ Ut = -K.[/. - yux + §([/, [/)[/., 

\ut = Uxxx - lu^U^ + \u^{U, U) + M(f/, f/^) + ([/, C4^) + ([/^, Ux), 

\ Ut = -K.[/. - \u^Ux + c/)c/. + {U, Ux)U, 

\ Ut = - \uSix + \u..,{U. U) + tz(f/, f/,) + (f/, f/,,) + (?7,, Ux), 
\ Ut = UxxU + UxUx - uUxU - \u'^Ux + \{U, U)Ux + {U, Ux)U, 

( Ut = Uxxx - \u^Ux + \Ux^P, U) + ul^U, Ux) + (t/, t/xx) + (t^x, t^x) 

^W,u)\ 

I = - \u^Ux - \{U, U)Ux + \u{U, U)U, 

{Ut = Uxxx - \u^Ux + U) + U(C/, C/j;) + {U, Uxx) + (C^x, C^x) 

-y{U,U) + \{U,U)\ 
Ut = + \{U, Ux)U - \u^U + \u{U, U)U, 



(4.15) 
(4.16) 
(4.17) 
(4.18) 
(4.19) 

(4.20) 

(4.21) 

(4.22) 

(4.23) 

(4.24) 

(4.25) 
(4.26) 





Uxxx + UUxU + U^Ux + {U, U)Ux + (t/, 
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ut = Uxxx - Gu'^u^ + 6u,^{U, U) + 12u(f/, ^4), 



(4.27) 



Ut = u^xx - Gu^u^ + u^{U, U) + 2u{U, U^) + {U, U^^) + {U^, U^), 
Ut = U^xx + ^u^xU + 3u.Mx - QuuJJ - 3u^U^ + {U, U)U^ (4-28) 

Ut Uxxx 

Ut = Uxxx + GuxxU + 6uxUx - UuUxU - 6u^Ux + {U, U)Ux (4-29) 
+ A{U,Ux)U, 

Ut = Uxxx - Gu^Ux + Ux{U, U) + 2u{U, Ux) + (t/, Uxx) + {Ux, Ux), 

Ut = -2Uxxx - QuxxU - 6uxUx + l2uUxU + Qu'^Ux + (f/, U)Ux (4-30) 
-2{U,Ux)U. 

All systems (I4.6p - (j4.3(jp admit the reduction f/ = 0. From this viewpoint, (j4.6p - (j4.9j) . 
(ICTl - dTOll . (ITOIl - dCTll and (ICTD - (ICT|) are generalizations of the 3'"'^ order Burg- 
ers equation, the trivial equation ut = 0, the pKdV equation and the mKdV equation, 
respectively. Actually, ()4.6p - ()4.8|) are the 3'''^ order symmetries of the 2"^^ order systems 
(|0|) - (P3j) . respectively. 

4.2 Integrability of systems (I1131)-(|113n]) 
4.2.1 Systems (fOjl and (flTHjl 

We investigate the 2"'^ order system ()4.3|) and its 3'''^ order symmetry ()4.6|) together. We 
note that the linear term Uxx in ()4.3|) vanishes iff a = — |, while the term Uxxx in ()4.6|) 
vanishes iff a = 0. This indicates that the case distinctions of 07^ — |ora = — ^ and a 7^ 
or a = are not very essential for the whole hierarchy of systems starting from (j4.3|) . As 
an extension of the Hopf-Cole transformation, we consider the change of variables. 



yj = qI udx'^ 

1 rx , 



Then we can triangular linearize (|4.3|) and (j4.6|) simultaneously as 

Wt = \{l + 2a)wxx + l{W,W), 

Wt = Wxx. 

and 

Wt = awxxx + {W, W)x, 

Wt = Wxxx- 



(4.31) 



(4.32) 



For some values of a, we can solve these systems easily. When a = 1, we can fully 
linearize systems (|4.3ip and (j4.32p through defining new variables V and v by W = Vx, 
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w + ^{V,V) = V (see [17]). When a = — |, we integrate the equation for w in (j4.3ip to 
obtain ^ 

w{x,t) =w{x,0) + - [ {W{x,t'),W{x,t'))dt'. 
3 Jo 

Similarly, when a = 0, we can integrate the equation for w in ()4.32|1 . We mention that 
Beukers, Sanders and Wang [37, 38] studied higher symmetries of the triangular linear 
systems fl4.31|l and ()4.32|1 in the case of scalar W. 

Symmetrization. We discuss symmetrization for the 2"^^ order system ()4.3j) . since it is more 
fundamental than its 3'''^ order symmetry ()4.6|) . To identify ()4.3p as a multi-component 
generalization of a system in [18,26], we assume the condition a ^ — ^ and rescale variables 
as 

dt = \{l + 2a)ds, u = Au, U = VQU'. 
In addition, we introduce a new parameter a by the relation 

3 

1 - 2a, 



1 + 2a 

where a 7^ |. Then fl4.3|) is rewritten as 

< = <. + 8w'< + (2-4a)(f/',f/i), 

U', = {1- 2a)U'^^ - Aau'^U' + (4 - 8a)u'U'^ - (4 + 8a)u'^U' 
- {2-Aa){U\U')U'. 

In the case where U' is scalar, this system is identical to (3.7) in [18]. In that case, 
considering the linear change of variables 

u' = q + r, U' = q — r, 

we obtain a system of two symmetrically coupled Burgers equations, which coincides with 
(3.6) in [18] or (3.7) in [26]. We note that system ()4.3|1 with a = — | does not have any 
counterpart in [18,26], because of its degeneracy of the linear part (cf. the introduction). 

4.2.2 Systems dHH) and P?f|> 

We investigate the 2"^^ order system ()4.4|) and its 3'''^ order symmetry ()4.7|) together. Here, 
we note that ()4.4|1 and ()4.7|1 are obtained from ()4.3|1 and ()4.6|1 . respectively, by rescaling t, 
U appropriately and taking the limit a — 00. Then, through the same change of variables 
as in section E. 2. H 

we can transform systems ()4.4|1 and ()4.7|1 to 

Wt = w^, + {W,W), 
Wt = 0, 
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and 

wt = w^xx + {W, W) 



(4.34) 

Wt = 0. 

Moreover, introducing a function g{x) such that g"{x) = {W,W), we can hnearize the 
equations for w in ()4.33p and ()4.34|) with respect to the variable w + g{x). 

Symmetrization. We discuss symmetrization for the 2"^^ order system ()4.4|) . To identify 
()4.4|) as a muhi- component generahzation of a system in [18,26], we rescale the dependent 
variables as 

u = 4m', U = 2V5U'. 

Then ()4.4|) is rewritten as 

u', = u'^^ + 8u'u'^ + 10{U',U'J, 
Ul = -2<f/' - Su'^U' - 10(f/', U')U'. 

In the case where U' is scalar, this system should coincide with (3.10) in [18], if it were 
written correctly. Unfortunately, in [18], Foursov made a mistake in deriving the equation 
for z in (3.10) from (3.9). It should be corrected as Zt = —2wxZ — Sw'^z — lOz^. If we 
consider the linear change of variables 

u' = q + r, U' = q — r, 

we obtain a system of two symmetrically coupled Burgers equations, which coincides with 
(3.9) in [18] or (3.6) in [26]. 

4.2.3 Systems dUSD and (HTHj) 

We concentrate our attention on the 2^^*^ order system ()4.5|) . and do not study its 3'''^ order 
symmetry ()4.8|) . Defining a new variable w hj w = ^{U, U), we find that ()4.5|) contains a 
two-component Burgers system, 

Ut = U^x + 2uUx + Wx-, , ^ 

, , (4.35) 

Wt = [uw)x. 

Therefore, system (14. 5p is a triangular system that consists of the Burgers system (j4.35|) 
and the linear equation for U with Burgers-system-dependent coefficients. We mention 
that in the long- wave limit (disappearance of Uxx), ()4.35p reduces to the Leroux system 
and that ()4.35p can be rewritten as a non-evolutionary scalar equation in (at least) two 
different ways. The symmetry integrability of (j4.35p as well as the existence of a recursion 
operator has already been demonstrated [18,62]. Nevertheless, we could find neither a 
linearizing transformation nor a true Lax representation for ()4.35p . In what follows, we 
discuss travelling- wave solutions of ()4.35|) . which are expected to give useful information 
on its properties. 

Substituting the travelling-wave form 

u{x, t) = f{z) — a, w{x, t) = g{z), z = x — at, 
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into (j4.35p . we get a system of two ordinary differential equations. Integrating it once, we 
obtain 



f' + p-af + g + b = 0, 
fg + c = 0. 



(4.36) 



Here, b and c are integration constants that are determined from the boundary conditions 
for u and w. Plunging g = —c/f into the first equation in ()4.36|) . we obtain the ordinary 
differential equation for /, 

d/ ^ _f-aP + V-c^ 
dz f 

For the sake of simplicity, we assume that — a/^ + bf — c can be factorized into the 
product (/ — ai){f — a2)(/ — c^s) with three distinct real roots ai, 02, 03. Thus we have 

a = «! + 02 + Q^3; b = a\a2 + «2a3 + a3ai, c = aio;2«3- 

Furthermore, we assume the conditions 7^ (j = 1,2,3) to obtain nontrivial solutions 
of ()4.35p . Indeed, if = 0, then c = and we obtain from ()4.36|) either a trivial solution 
or a solution of the scalar Burgers equation. Noting the identity 



/ 



(/ - "i)(/ - "2)(/ - «3) ("1 - a2)(ai - «3) V/ - "1 / - "3 



+ 



02 



[a2 - ai){a2 - 03) V / - "2 / - as 



we can integrate ()4.37|1 to obtain 



1 + 



tt3 — Ctl \ ("l-«2)(«l-a3) 

/ - "3 



1 + 



^3 - "2 
/ - "3 



(a!2-"l){"2-°=3) 



de 



(4.38) 



where c? is a constant. Now, it is clear that the functional form of !/(/ — 0:3) depends on 
the ratio of two powers on the left-hand side, (og^ — a^^)/(a]~^ — 0:3^). Let us consider 
the simplest case in which this ratio is unity, i.e.. 



1 

03 



1 1 
2 \ai a2 



In this case, we have ai + 02 7^ for the existence of as. Then we can solve ()4.38j) 
explicitly for !/(/ — a^): 



1 



«i + «2 + ^^A/l + exp 



f -as 

Here, Zq is the constant given by 



2a;ia;2 



(4.39) 



— 4aia2 
(«! + 0:2)^ 



(01-12) 
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Using the relation = 2aia2/{ai + 02), we can rewrite (j4.39|) as 



m = ^ ^ . (4.40) 



2 



In order for the function f{z) to be non-singular, we should assume the condition ai{ai + 
^2) > 0. To summarize, we have obtained in the simplest case a travelling-wave solution 
of ()4.35p given by 

u{x, t) = f{x — at) — a, w{x, t) 



fix-aty 

with ()4.4()j) . a = ai + a2 + 2aia2/ (ai + 02) and c = 2(aia2)^/ (^i + «2)- 

When a nontrivial solution, like the above, of subsystem ()4.35|) is given, we can solve 
the remaining equation for U in the original system ()4.5|) . {Uj)^ = (uUj)^, in the same 
way as in section IB. 2. 31 

Symmetrization. We discuss symmetrization for the 2^^*^ order system ()4.5|) . With the 
rescaling of dependent variables 

u = 2u', U = \^U', 

()4.5|1 is rewritten as 

u[ = u',, + 4u'u', + 3{U',U'J, 
U[ = u'JJ' + 2u'U'^. 

In the case where U' is scalar, this system is identical to (3.4) in [18]. If we consider the 
linear change of variables 

■u' = g + r, U' = q — r, 

we obtain a system of two symmetrically coupled Burgers equations, which coincides with 
(3.3) in [18] or (3.5) in [26]. 



4.2.4 System (glU) 

In the case where U is scalar, system (j4.9j) coincides with system (|4.8p . However, unlike 
()4.8|) , ()4.9p in the general N case does not possess a 2"^^ order symmetry of the form 
()4.1|) . System ()4.9|) contains the 3'''^ order symmetry of the two-component Burgers sys- 
tem ()4.35|) . where w is again given by w = ^{U,U). Therefore, in order to demonstrate 
the integrability of (|4.9p . we first of all need to know either a linearizing transformation 
or a Lax representation for (j4.35j) . This remains as an open question. 

Symmetrization. Through symmetrization of (|4.9p in the case of scalar U, we just ob- 
tain the 3'''^ order symmetry of the two symmetrically coupled Burgers equations in sec- 
tionlOSl 
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4.2.5 System 

System ()4.10|) is connected with system ()4.12|) through a Miura-type transformation. We 
discuss the integrabihty of these two systems in section 14. 2. 71 

Symmetrization. In the case where U is scalar, we consider the hnear change of variables 

u = q + r, U = \/a{q — r). 

Here, a is a nonzero constant. Then we can rewrite ()4.10|) as a system of two symmetrically 
coupled equations, 

qt = \qxxx - \r^^^ + 1(1 + 3a) (g - r)g,, + |(1 - 3a) (g - r)r,, 
+ hi - + 3a(g - rfr.^ - |a2(g - r)^ 

'^t = -\qxxx + \rxxx - 1(1 - 3a)(g - r)^,, - \{l + 3a)(g - r)r^^ 
- \ql + \rl + 3a(g - r^q^ - \a^{q - rY- 

This system does not belong to the class of systems studied in [25,26], because of its 
degeneracy of the linear part. 

4.2.6 System dUTTD 

For system ()4.11|) . we have the relation (m, — (f/, U))t = 0. Thus, we can set 

u,-{U,U) = <j){x), (4.41) 

where the function (f){x) does not depend on t. Then the equation for U is rewritten in 
terms of (j){x) as 

Ut = U,,, + 2(f)U, + (j),U. (4.42) 
The solutions of ()4.42|) are given by 

where \l/(x; A) is a solution of the ordinary differential equation 

+ 20^, + 0,^ = A^. (4.43) 

Once we obtain 0(x) and U{x,t), we can determine u{x,t) by using 1)4.4111 . The vec- 
tor equation ()4.43j) is of the same form as the scattering problem associated with the 
Kaup-Kupershmidt equation [63,64]. We can see that this is not a coincidence through 
investigation of the 5^^ order symmetry of system 1)4.111) . Indeed, the 5^^ order symmetry 
is rewritten (up to a scaling of t^) in terms of and \E' as 

0*5 + 4>xxxxx + 1O00XXX + 250,0,, + 200^0, = 0, 

= 9^,,,,, + 300^^,,, + 450,*,, + (350,, + 2002)*, 
+ (100,,, + 2000,)*. 
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The first equation is nothing but the Kaup-Kupershmidt equation, while the second 
equation together with ()4.43|) constitutes a Lax representation for it. 

Symmetrization. In the case where U is scalar, we consider the linear change of variables 

u = q + r, U = y/a{q — r). 

Here, a is a nonzero constant. Then we can rewrite ()4.11|) as a system of two symmetrically 
coupled equations. 



qt = \qxxx - \rxxx + {\ + oi){q- r)q^^ + {\-a){q- r)r^x 
+ (1 - \a)ql + aq^r^ - (1 + \a)rl - a{q - rfq^ 
+ 3«(g — tYtx — oP'iq — t)^, 

n = -\qxxx + - (I - a) (g - r)q^^ - {\ + a){q- r)r^^ 

- (1 + \a)ql + aq^r^ + (1 - \oi)rl + 3a(g - rfq^ 

— a{q — r)^r^ — a'^{q — r)"^. 

This system does not belong to the class of nondegenerate systems studied in [25,26]. 
4.2.7 System dHHD 

For system ()4.12|) . if we define new variables w and W by 

W = U^ + uU, 

they satisfy the following system: 

wt = -3{W,W^), 
Wt = W^^^ + w^W + 2wW^. 

This system coincides with the multi-component Drinfel'd-Sokolov system ()3.3p . up to a 
scaling of W. The Miura map ()4.44p is a multi-component generalization of that for the 
case of scalar U in [43,44]. 



Relation to system { 4-10^ - If we introduce a new scalar variable v by 

V = u^~u^ + {U,U), (4.46) 



system (I4.12p is changed into the following system: 

vt = {3v{U, U) + 3(f/, U^^) - 3{U, Uf)^, 
Ut = a^xx + v^U + vU^ - 3{U, U,)U. 



(4.47) 



Then, it is straightforward to obtain ()4.10|) (for u and U) from ()4.47p through potenti- 
ation V = Ux- Combining (j4.46j) and (j4.44p . we obtain the relation v — 2w = 3Ux, and 
consequently, 

PX 

M — 2 / w dx' = 3m. 
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Using this relation, we can also rewrite (j4.44|) as a transformation between system (j4.1L)j) 
and the multi-component Drinfel'd-Sokolov system ()4.45|) . 



Symmetrization. In the case where U is scalar, we consider the linear change of variables 

u = q + r, U = ^/oi{q — r). 

Here, a is a nonzero constant. Then we can rewrite ()4.12|1 as a system of two symmetrically 
coupled equations. 



qt = \[qxx - rxx + (1 + - r)q^ + (1 - - r)r^ 

- (1 - _ (1 + + (1 - a)qr^ + (1 + ayi, 

(4.4oj 



rt = + - (1 - a) (g - r)g^ - (1 + a)(g - r)r^ 

+ (1 + a)q^ + (1 - a)q'^r - (1 + a)qr'^ - (1 - Ci)r^]x- 

This system does not belong to the class of nondegenerate systems studied in [25,26]. 
However, it was found in connection with the Kac-Moody Lie algebras and written in a 
Hamiltonian form about twenty years ago [43,44]. More specifically, system ()4.48p with 
a = — 1 coincides with (3)-(4) in [43] for the Dg^-* case, up to a scaling of variables (see 
also the generalized mKdV equation in [44] for the Ag^'* case). 

4.2.8 System dHHD 

System (|4.13p is merely a potential form of the multi-component Zakharov-Ito system 

Symmetrization. In the case where U is scalar, we set 

u = q + r, U = q — r. 
Then we can rewrite fl4.13|) as a system of two symmetrically coupled pKdV equations, 

Qt 2Q.XXX 2^ XXX ~l~ '^Qx ~^ ^x) 

2Q.XXX 2^xxx ~^ Qx ~^ '^^x' 

This system is identical to (37) in [25] or (3.10) in [26]. 

4.2.9 System dHHD 

Remark. If we consider separately the two cases a = and a 7^ 0, we can reduce the num- 
ber of parameters in system ()4.14j) by scaling variables. In the former case the parameter 
h can also be fixed at any nonzero value, while in the latter case only the parameter a 
can be scaled away. However, this case distinction is neither necessary nor essential, as is 
demonstrated below. 

If we define new variables w and W by 

w = u^ + ^{U,U), 
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they satisfy the following system: 



wt = w^^^ + Qww^ + - t) ^)^' (A 4g^ 

Wt = 2{wW)^. 

Thus the essential parameter is 6 — rather than a or b. If 6 — 7^ 0, system 
fl4.49|) coincides with the multi-component Zakharov-Ito system ()3.5|) . up to a scaling 
of variables. When b — a^/4 = 0, ()4.49|) is a triangular system that consists of the 
KdV equation and the linear equation for W with KdV-equation-dependent coefficients. 
This triangular system was studied from a point of view of symmetries in [65] (see also 
[22,66,67]). As we have noted in section I3.2.3t we can relate with W the inverse square 
of a solution of the KdV linear problem. 

Symmetrization. In the case where U is scalar, we set 

u = q + r, U = q — r, a = 1 + 2a, b = 2/3. 

Then we can rewrite (I4.14|) as a system of two symmetrically coupled pKdV equations. 



qt = \qxxx + \r^xx + (1 + a)(g - r)q^^ - a{q - r)r^^ + (3 + a)ql 
+ (2 - 2a)q^r^ + (1 + a)rl + (2 + 4a) (g - r^q^ + p{q - r)^ 

n = hxxx + hrxxx + Oi{q - r)qxx - {I + a){q - r)r^^ + (1 + a)ql 



(4.50) 



+ (2 - 2a)qxrx + (3 + a)rl + (2 + 4a) (g - r fr^ + [3{q 



r 



i4 



This system is identical to (12) in [25]. Foursov claimed in [25] that this system was 
either reduced to the representative case a = (a = 1) or decoupled by a linear change 
of dependent variables. However, in fact this is not true. As far as we consider a linear 
change of variables, we need one more representative case, a = — | (a = 0), in which 
(j4.5(J|) cannot be decoupled. Therefore, we can say that (at least) one system is missing 
from the final list of Foursov-Olver given in [26]. 

4.2.10 System dHHD 

System ()4.15|1 is merely a potential form of the Jordan KdV system ()3.4j) . We see in 
section l4. 2. 221 that this system is also connected with system ()4.27|1 through a Miura-type 
transformation. 

Symmetrization. In the case where U is scalar, we consider the linear change of variables 

u = ^(g + r), U = ]^{q-r). 

Then ()4.15p is decoupled into two pKdV equations, 

qt = qxxx + 3g^, 

I t — ' XXX "T Ol ^. 

This corresponds to (27) in [25]. 
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4.2.11 System (jUISD 

System ()4.16|) is connected with system ()4.28p through a Miura-type transformation. We 
discuss the integrabihty of these two systems in section 14.2.231 

Symmetrization. In the case where U is scalar, we consider the hnear change of variables 

u = q + r, U = V3{q — r). 

Then we can rewrite (I4.16p as a system of two symmetrically coupled pKdV equations, 

qt = Qxxx + 3(g - r)q^^ + 3{q^ + r^)q^ + 3(g - rfq^, 
rt = r^^^ - 3(g - r)r^^ + 3(g^. + r^)r^ + 3(g - r)V^. 

This system coincides with (34) in [25] or (3.9) in [26]. 

4.2.12 System (flTTj) 

System ()4.17|) is connected with system ()4.29p through a Miura-type transformation. We 
discuss the integrabihty of these two systems in section 14.2.241 

Symmetrization. In the case where U is scalar, we consider the linear change of variables 

u = -{q + r), U = —{q-r). 

Then we can rewrite (I4.17|) as a system of two symmetrically coupled pKdV equations, 

Qt = Qxxx + 3(g - r)q^^ + 3ql + 3(g - rfq^, 
n = r^xx - 3(g - r)r^^ + 3rl + 3{q - r fr^. 

This system coincides with (28) in [25] or (3.8) in [26]. 

4.2.13 System (HHED 

System ()4.18|) is connected with system ()4.30p through a Miura-type transformation. We 
discuss the integrabihty of these two systems in section 14.2.251 

Symmetrization. In the case where U is scalar, we consider the linear change of variables 

u = q + r, U = V3{q — r). 

Then we can rewrite (I4.18p as a system of two symmetrically coupled pKdV equations, 

Qt = -\qxxx + frxxx - 6(g - r)r^^ + Qrl + 12(g - rfr^ 
+ 3(g-r)^ 

rt = \qxxx - \rxxx + 6(g - r)q^^ + 6ql + 12{q - rfq^ 
+ 3{q-r)\ 

This system is identical to (13) in [25] with a = 6. Thus, it is equivalent to (41) in [25] 
or (3.12) in [26], up to a linear change of dependent variables. 
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4.2.14 System (jUinD 

For system ()4.19p . if we introduce a new variable w by 

w = u^ + 2{U,U), (4.51) 

it solves the KdV equation, 

Wt = Wxxx + '2.wWx- (4.52a) 

Therefore, system ()4.19p is reduced to a triangular form. The equation for U is rewritten 
in terms of w as 

Ut = 4U^^^ + w^U + 2wU^. (4.52b) 

We note that the vector equation ()4.52b|) is of the same form as the time part of the 
linear problem for the KdV equation (j4.52a|) . This relation between system (|4.19p and the 
KdV equation resembles the relation between the 5*^ order symmetry of system (|4.1ip and 
the Kaup-Kupershmidt equation shown in section 14.2.61 A recursion operator and a Lax 
representation for the triangular system ()4.52|) were given in [67] and [24], respectively. 
Once we obtain w{x,t) and U{x,t), we can determine u{x,t) by using 1)4.511) . 

Symmetrization. In the case where U is scalar, we consider the linear change of variables 

u = 6{q + r), U = V3i{q — r). 

Then we can rewrite (I4.19|l as a system of two symmetrically coupled pKdV equations, 

qt = \qxxx - \rxxx + 6(g - r)ga;a; + 6g^ + Ylq^^^r^ - 6rl 
-12(g-r)2(g,-r,)-3(g-r)^ 

n = -\qxxx + I'^xxx - 6(g - r)r^^ - Qql + 12g^r^. + 6r^ 
+ 12(g-r)2(g,-r,)-3(g-r)^ 

This system coincides with (42) in [25] or (3.13) in [26]. 

4.2.15 System (gHOD 

For system ()4.2flj) . if we introduce a new variable w by 

w = -u^ + ]^u'-]^{U,U), (4.53) 

it solves the KdV equation, 

= Wxxx — 'iwWx- (4.54a) 

Therefore, system ()4.2U|) is reduced to a triangular form. Still, it is a very interesting 
system. To see this, we rewrite the equation for u in terms of w as 

Ut = —^ul + v^Ux — "iwUx — — Wxx- (4.54b) 

Then, the reduction w = changes this equation to a nontrivial closed equation for u. 
With a rescaling of variables, it reads 

Ut = Ux{ux-u^). (4.55) 
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Equation (j4.55p possesses an infinite set of commuting symmetries 

Mr„ = Ux{ux - u^Y, n G M. 

We can easily obtain a travelling- wave solution of (j4.55j) with two arbitrary constants [68], 
which is called a complete solution in the theory of partial differential equations. However, 
we do not know any explicit formula for the general solution of fj4.55|) . Using ()4.53p . we 
can rewrite the equation for [/ as a linear equation with a m, w-dependent coefficient. 

Symmetrization. In the case where U is scalar, we consider the linear change of variables 

u = q + r, U = q — r. 

Then we can rewrite fl4.2U|) as a system of two symmetrically coupled mKdV equations, 

Qt = \(lxxx + \^xxx ~ ^){(lxx ~ ^xx) ~ {Qx ~ ^x)fx 

^ +{q^ - 5qr)qx - {q^ + qr)rx, 

2Q.XXX 'qJ'xxx ~^ ^^(^qxx ^xx^ ~l~ (Q'i' ^x^qx 

— {qr + r'^)qx + (— 5gr + r'^)rx. 
This system coincides with (57) in [25] or (3.18) in [26], up to a scaling of variables. 

4.2.16 System (|1I2H) 

In system (14.211) . u and {U, U) satisfy the same equations as in system ()4.20p . Thus, if we 
define w by ()4.53p . we obtain ()4.54a|) and ()4.54b|) again. The only difference between the 
two systems ()4.21|1 and ()4.2()j) lies in the forms of equations for f/, which can be rewritten 
as linear equations with u, t(7-dependent coefficients. 

Symmetrization. In the case where U is scalar, 1)4.211) is identical to ()4.2()|1 . Thus, through 
symmetrization, we just obtain the same result as in section 14.2.151 

4.2.17 System 

System ()4.22|) has already been obtained in [69] as a reduction of a bi-Hamiltonian system 
(see also [70]). If we introduce a new variable w by 

w = -Ux + - ^(f/,t/), 

it solves the KdV equation, 

Wt = Wxxx — SwWx- (4.56) 

Therefore, system ()4.22|1 is reduced to a triangular form. Substituting \{U,U) = —u,j. + 
— w into the equations for u and U respectively, we obtain 

r Mt = -{WU + Wx)x, 

\ut = -{wU)x. 
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Thus the reduction = is trivial in this system. We mention again (cf. section |4.2.9|) 
that the above equation for U coupled to the KdV equation ()4.56|) was studied in [65]. 



Remark. Actually, a one-parameter deformation of system ()4.22j) . 

(4.57) 



Ut = 3aux + u^^^ - \u^u^ + \u^l^, U) + u{U, U^) + (f/, U^^) + f4), 



Ut = aU^ + u^^U + u^U^ - uu^U - ^u^U^ + |(f/, U)U^ + {U, U^)U, 

is still integrable. Indeed, if we introduce w in this case by 

1 1 

w = -a-u^ + -u^- -{U, U), 

system fl4.57|) is changed into the following system: 

wt = uixxx - Sww^ + 2a{U, U^), 
Ut = -{wU),. 

When a 7^ 0, this system coincides with the multi-component Zakharov-Ito system ()3.5p . 
up to a scaling of variables. 

Symmetrization. In the case where U is scalar, we consider the linear change of variables 

u = q + r, U = q — r. 
Then we can rewrite fl4.22|) as a system of two symmetrically coupled mKdV equations, 

Qt = \qxxx + \rxxx + {q- r)qxx + {Qx - rx)qx - ^rq^ - 2q'^r^, 
n = \qxxx + \r^^^ - (g - r)r^^ - {q^ - r^)r^ - 2r^q^ - Aqrr^. 

This system coincides with (58) in [25] or (3.19) in [26]. 
4.2.18 System (jOSD 

For system ()4.23|) . if we introduce a new variable w by 

w = -u^ + - ^{U,U), 

it satisfies the KdV equation, 

Wt = Wxxx — SwWx- (4.58a) 

Therefore, system fl4.23|) is reduced to a triangular form. The equation for u is rewritten 
in terms of w as 

«,^-A..i„'-... + »..-...-2.«= + 2.=. (4.58b) 

The triangular system ()4.58|) possesses (at least) two higher symmetries. The first higher 
symmetry is given by 

Wt^ = Wxxxxx - 5ww^^^ - lOw^w^^ + yU'^Wx, 

Ut^ = —\u'^W + 2u'^w'^ — 2w^ + U^WUx — \w'^Ux — U^Wx + 5uWWx 
+ 2uUxWx + Swl — U^Wxx + 5wWxx + UxW 

XX XXX ^ XXXX^i 
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which obviously vanishes under the reduction w = 0. Similarly, the second one vanishes 
under the same reduction. On the other hand, the reduction w = changes ()4.58p to a 
nontrivial closed equation for u, 

2 I 4 
Ut = —U Ux + -u . 

As far as we could check with the help of a computer, this equation seems to have no 
polynomial higher symmetry. We can construct its general solution in implicit form using 
the method of characteristic curves. 

Symmetrization. In the case where U is scalar, we consider the linear change of variables 

M = g + r, U = q — r. 
Then we can rewrite (I4.23j) as a system of two symmetrically coupled mKdV equations, 

Qt '2QXXX ~l~ 2^xxx ~l~ ^)(^Qxx ^xx) i^Qx ^x^^x 

— (3gr + r'^)qx + (— 3gr + r'^)rx + ^{q — r)^q, 

ft 2Q.XXX ~l~ 2'^xxx ~l~ f)(^Qxx ^xx) ~l~ i^Qx ^x^Qx 

+ {q"^ ~ 3gr)ga; — (g^ + 3qr)rx — ^{q — r)^r. 
This system coincides with (61) in [25] or (3.20) in [26], up to a scaling of variables. 

4.2.19 System dOH) 

For system ()4.24j) . if we introduce a new variable w by 

w = -u, + ^u^-^{U,U), (4.59) 

it satisfies the KdV equation, 

Wt = Wxxx — SwWx- 

Therefore, system ()4.24|1 is rewritten in a triangular form. Substituting |(f/, U) = —u^ + 
— w into the equations for u and U respectively, we obtain 

Ut = —WxU — }:WV? — Wxx + w'^ , 

Ut = -2(w^x + wu)U. 

Thus the reduction w = is trivial in this system. We note that in ()4.60|) no term 
involves x-derivatives of u, U such as Ux, Ux, Uxx, Uxx- Then, for a given solution of the 
KdV equation w(x, t), the equation for m(x, t) can be regarded as a Riccati equation with 
X fixed. Once we obtain w{x,t) and u{x,t), we can integrate the equation for U as 

U{x,t) =e-5/o("'-+"'")'i*'f/(x,0). 

Remark. Actually, system ()4.24|) is the 3'''^ order symmetry of a nontrivial 1*^* order system, 

Ut,=Ux + l{U,U), 
Ut, = \uU. 
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For this system, w defined by (I4.59p obeys the hnear equation Wt^ = and the equation 
for u is rewritten as Ut^ = ^u^ — w. System ()4.61|) in the = 1 well as its higher 

symmetries was studied in [56,71,72] (see also [73] for its soliton-like solutions). 

Symmetrization. In the case where U is scalar, we consider the linear change of variables 

u = q + r, U = q — r. 

Then we can rewrite (I4.24j) as a system of two symmetrically coupled mKdV equations, 

qt = 2^xxx "I" 2'^xxx "I" 4(0' 1^)i'i(lxx "^xx) "I" 2(^3; "^x)"^ 

+ — 6qr — r^)qx — {q^ + 2qr)rx — q^r{q — r), 

rt = \qxxx + \rxxx + \ {q - r){qxx - 3r^^) + - r^Y 
— {2qr + r'^)qx — ^{q"^ + Qqr — r^)r^ + qr^{q — r). 

This system coincides with (23) in [25] with a = 1, up to a scaling of dependent variables. 
Thus, it is equivalent to (62) in [25] or (3.21) in [26], up to a linear change of variables. 

4.2.20 System (glSlD 

System ()4.25p is just a disguised form of a single vector equation. Indeed, if we introduce 
an (A^ + l)-component vector W hj W = {u, U) = {u, Ui, . . . , U^), system (j4.25j) can be 
rewritten in the form 

Wt = Wxxx + {W,W)Wx. 

This is a well-known vector mKdV equation and its integrability has been established in 
the literature [27,28,74,75]. 

Symmetrization. In the case where U is scalar, we consider the linear change of variables 

u = i(g + r), U = ^{q-r). 

Then we can rewrite ()4.25p as a system of two symmetrically coupled mKdV equations, 

j qt = Qxxx + qrqx, 
\rt = Txxx + qrvx- 

This system is known as (the non-reduced form of) the complex mKdV equation [76]. It 
is identical to (43) in [25] or, (3.14) in [26] with a correction of misprints. 

4.2.21 System 

System ()4.26|1 is also a disguised form of a single vector equation. Indeed, if we introduce 
an (A^ + l)-component vector W hy W = {u, U) = {u, Ui, . . . , Un), system ()4.26j) can be 
rewritten in the form 

Wt = Wxxx + (W, W)Wx + {W, Wx)W. 
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This is another well-known vector mKdV equation, for which a Lax representation was 
given in [77] for the = 1 case and in [78] for the general case. 

Symmetrization. In the case where U is scalar, we consider the linear change of variables 

u = q + r, U = i{q — r). 
Then we can rewrite (I4.26p as a system of two symmetrically coupled mKdV equations, 

qt = q-xxx + Qqrqx + '^q'^r^, 
n = r^^^ + 2r^q^ + 6qrr^. 
This system coincides with (48) in [25] or (3.15) in [26], up to a scaling of variables. 

4.2.22 System dHHD 

System ()4.27|1 is known as a Jordan mKdV system [27]. Let us briefly summarize its 
integrability. It is well-known that the matrix mKdV equation, 

Qt = Qxxx - 3(g,Q2 ^ g2g^) ^ ^4_g3) 

admits a Lax representation [52,75,78,79]. Then, system ()4.27|) is also integrable, because 
it is obtained from ()4.63j) through the following reduction: 

N 

Q = ul + UjBj, {e-i, ej}+ = -2Sijl. 
i=i 

We mention that ()4.27p admits a generalization to a system for two vector unknowns 
preserving the integrability [29]. 



Relation to systems (^^) and { 4-l^ . If we define new variables w and W by [27] 

w = ±u^-u^ + {U,U), 



W = U^T 2uU, 
they satisfy the following system: 

wt = w,^^ + 3{w^ -{W,W))^, 

This system coincides with the Jordan KdV system ()3.4|) and, through potentiation of it, 
we obtain system (j4.15|) . 

Symmetrization. In the case where U is scalar, we consider the linear change of variables 

M = i(g + r), U = ^{q-r). 
Then ()4.27p is decoupled into two mKdV equations, 

qt = qxxx - Qq'^qx, 



' t ' XXX ^' ' X- 



This corresponds to (50) in [25]. 
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4.2.23 System (jOHD 

We note that through introduction of a new scalar variable w by 

,2 



W =Un 



u 



system fl4.28|) is transformed to a system of coupled KdV-mKdV type, 



(4.64) 



Let us demonstrate that system ()4.28|) admits a Lax representation. We consider a pair 
of linear equations for a column- vector function -0, 

with the matrices U and V of the following form: 



U 



( -4;iCHi - 2iCQR 

+ Q^R - QR, + 2QPR 



V 



-Kh Q 
R Kim + P 



ACQ + 2iC(Q. + QP) - Qxx 
- 2Q,P - QP, + 2QRQ - QP^ 



(4.65a) 



4CR + 2iC{-Rx + PR) 

— Rxx + PxR + 2PRx 

+ 2RQR - P'^R 



\ 



(4.65b) 



4iC'/™ + 2iCRQ - Pxx + RxQ 
- RQx + P^P - PP, + 2PRQ 
\ +2RQR- P'^R +2RQP - P^ -Sg^P J 

Here, ( is the spectral parameter, J/ and Im are the / x / and m x m unit matrices 
respectively, Q, R and P are I x m, m x I and m x m matrices respectively, and g 
is a scalar function. The compatibility condition -0^^ = i/?^^ implies the so-called zero- 
curvature condition, 

Ut-v^ + uv -vu = o. 

Then, substituting ()4.65|) into this condition, we obtain a system of three matrix equations, 

Qt + Qxxx + ^{QxP)x — ^QxRQ — ^QRQx + ^QxP^ 

+ 3gp,.p - 3(?,gp = o, 

-Rt + Rxxx — 'i{PRx)x — SRxQR — 3RQRx + SP'^R^ 
+ 3PPxR + 3gxPR = O, 

Pt + Pxxx + K9^P)x - SiPRQ + RQP)x + ^PPxP 
+ W^RQ - 3RQP'^ = O. 



(4.66) 
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We note that this system admits the reduction R = *Q, *P = —P, where the superscript 
* denotes the transposition. In particular, if we choose 

(Q = {u,W^,...,WN) = {u,W), 
/ u \ 



R 



P 



\Wn J 
( 



u 



Vi ■■■ \ 

o 



\ -Vn 

system (I4.66p is reduced to the system 



V 

-W o 



J 



(4.67) 



f ut + u^,, - Qu'^u^ - -iu^m, W) + {V, V)) - 3m((H^, W,) + (V, K)) 
+ ?,g,{W,V)-^W,,V), = Q, 

-?>{W,W,)W -?>{W,V),V = 

Vt + + Ka.v), - Ku'V), - 3((H^, V)W), - 3(V, V^)V 
- ?>u{V, V)W + 3u(W, V)V = 0, 

together with one constraint, 

{u 'VW - u 'WV)^ - {W, V) Cvw - 'WV) = O. 

When W and V are scalar, i.e. = 1, this constraint is satisfied automatically and we 
obtain a three-component mKdV system. There is another case, the case W = V, for 
which the constraint is satisfied. Then, if we set 



g = U, W = V= -^U, 

and change the sign of time t {t ^ —t), system ()4.(j7p collapses to system ()4.28j) . 
Relation to system { 4- 161 . If "we introduce a new scalar variable v by 

V = u^-u^ + ^{U,U), 

system fl4.28|l is changed into the following system: 

Vt = {v^x + Sv^ + v{U, U) + {U, f/^^))^, 
Ut = f/.xx + ^v^U + 3t;t/, + (f/, U^)U. 



(4.68) 



Then, it is straightforward to obtain ()4.16|) (for u and U) from ()4.68|) through potentiation 
V = Ux- It should be mentioned here that the authors encountered two papers [80,81] on 
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the integrability of (j4.68j) in the case of scalar U, after they had obtained all the presented 
results independently. 

Symmetrization. In the case where U is scalar, we consider the linear change of variables 

u = q + r, U = V3{q — r). 

Then we can rewrite fl4.28|) as a system of two symmetrically coupled mKdV equations, 

qt = [qxx + 3(g - r)q^ + q^ - 12gV + Sgr^]^, 
n = [txx - 3(g - r)r^ + Sg^r - 12qr'^ + r^]^. 

This system coincides with (55) in [25] or (3.17) in [26]. Moreover, if we introduce new 
variables q and f by 

q = VSige^"^'?-^)'^^', f = VSire-^'^^''-")'^^', 
they satisfy the system of coupled mKdV equations ()4.fj2j) . 

4.2.24 System (IQ^ 

Through introduction of a new scalar variable w by 

W = Ux — M^, 

system fl4.29|) is transformed to a system that looks very similar to ()4.64|) . 



Ut = U^xx + Q{wU), + (f/, U)U^ + 2{U, U),U. 



(4.69) 



System ()4.69p is a mult i- component generalization of a flow of the Jaulent-Miodek hier- 
archy [82]. Let us demonstrate that ()4.69|) admits a Lax representation. We consider a 
pair of linear equations for a column- vector function '0, 

i^,^ + {Q + Ci?) V = cv, 

i^, = {ACI + 2CR + 2Q + |i?2)V;, - [CR^ + Q^ + l{R%]i^. 

Here, ( is the spectral parameter, and Q and R are square matrices with the same dimen- 
sion. The compatibility condition -0^^^ = i/j^^^ for ()4.70|) implies a system of two matrix 
equations. 



Qt = Qxxx + 3(Q^)x + i{R'^)xxx + \R^Qx + \[Q{.R^)x + ^{,R^)xQ]-, 
Rt = Rxxx + S{QR + RQ)x + f [3(-R^)a;-R + R{R'^)x + '^R^Rx], 



together with one constraint. 



[Q,R'] = 0. 



39 



If we consider the reduction, 

N 

t). ' 

3~' 



Q = wl, R = — i ^ Ujej, {ei, ej}+ = -25ijl, 



the constraint is automatically satisfied and system ()4.71|) collapses to system ()4.69|) . This 
Lax representation for (j4.69p can be rewritten as that for (j4.29|) [83]. 



Relation to system ( -^. i? ). If we introduce a new scalar variable v by 

1 



v = u,-u^ + ^{U,U), (4.72) 





system (I4.29p is changed into the following system: 

vt = + Sv' + 2v{U, U) + (f/, U^,) + |(f/,, t/,.))^, 
Ut = U^xx + 6v.,U + 6vU, + 2{U, U,)U. 



(4.73) 



Then, it is straightforward to obtain ()4.17|1 (for u and U) from ()4.73|1 through potentiation 

V = Ux- 

Symmetrization. In the case where U is scalar, we consider the linear change of variables 

u = -{q + r), U = —[q-r). 

Then we can rewrite fl4.29|) as a system of two symmetrically coupled mKdV equations, 

Qt = [Qxx + 3(g - r)qx + - 6q^r + Sqr"^]^, 
ft = Vxx — "ii^q — + 3g^r — 6gr^ + r^]^. 

This system coincides with (51) in [25] or (3.16) in [26]. It is known as a fiow of the 
modified Jaulent-Miodek hierarchy [83] (see also (7.37) in [84]). While elaborating on 
this paper, the authors encountered one paper [85] on the three-component generalization 
of this fiow [ (irail with = 2]. 

4.2.25 System dUnOD 

For system ()4.30p . if we define new variables w and W by 

w^Ux + u' + \{U,U), 
W = Ux + 2uU, 

they satisfy the following system: 

Wt = Wxxx - 6wwx + 2{W, Wx), , ^ 

(4.75) 

Wt = -2Wxxx + QwWx- 

This system coincides with the multi-component Hirota-Satsuma system ()3.6|) . up to a 
scaling of variables. The Miura map (|4.74p is a multi-component generalization of that 
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for the case of scalar U in [43,44] and that for the case of two-component vector U in [61]. 



Relation to system d-^.i.^ . If we introduce a new scalar variable v by 



v = u^-u^-hu,U), (4.76) 
6 



system (I4.3(J|) is changed into the following system (cf. (4.3) in [86]): 

vt = {v.. + + Av{U, U) + 2(f/, U,^) + (t/,, U,) + |(f/, Uf)^, 
Ut = -2U^^^ - 6v^U - 6vU^ - 4(f/, U^)U. 



(4.77) 



Then, it is straightforward to obtain ()4.18j) (for u and U) from ()4.77|1 through potenti- 
ation V = Ux- Combining ()4.76p and ()4.74|) . we obtain the relation v + w = 2ux, and 
consequently, 

u + wdx' = 2u. 



Using this relation, we can also rewrite ()4.74|1 as a transformation between system ()4.18|1 
and the multi-component Hirota-Satsuma system ()4.75|1 . 

Symmetrization. In the case where U is scalar, we consider the linear change of variables 

u = --{q + r), U = —i{q-r). 
Then we can rewrite fl4.3()|l as a system of two symmetrically coupled mKdV equations, 

Qt = [-\qxx + \rxx + 3(g - r)g^. - 2r^]^, 
n = [hxx - \rxx - 3(g - r)r^ - 2(f'\^. 



This system is identical to (63) in [25] or (3.22) in [26]. It was found in connection with 
the Kac-Moody Lie algebras and written in a Hamiltonian form about twenty years ago 
(cf. the C^^^ case in [43] or the B^^^ case in [44]). 



41 



5 The case Xi = X2 = ^ — coupled Ibragimov— Shabat 
equations — 

In this section, we classify 2'^'^ order and order systems in the Ai = A2 = | (Ibragimov- 
Shabat weighting [36]) case. In the first part f section l5.1|) . we present a complete list of 
such systems with a specific order symmetry. In the second part (section l5.2|l . we prove 
that the listed systems are linearizable. 



5.1 List of systems with a higher symmetry 

The general ansatz for a Ai = A2 = | homogeneous evolutionary system of 2"^*^ order for 
a scalar function u and a vector function U takes the form 



= diUxx + ci2u'^Ux + 03^^ + a^UxiU, U) + a^u{U^ Ux) 

Ut^ = agUxx + aguujj + aioM^f4 + auu^U + ai2(f/, U)Ux 
+ ai3(f/, Ux)U + ai4M2(f/, U)U + ai^{U, UfU. 



(5.1) 



The following constraints guarantee the order to be 2 and the system not to be triangular: 

(ai, as) 7^ (0, 0), (04, 05, ae, 07) 7^ (0, 0, 0, 0), (ag, aio, an, 014) 7^ (0, 0, 0, 0). 
Similarly, the general ansatz for a 3'''^ order system takes the form 



(5.2) 



hiUxxx + h2U^Uxx + huul + hiU^Ux + hvJ + beUxxiU, U) 
+ hux^U, Ux) + hu{Ux, Ux) + hu{U, Uxx) + hou'^UxiU, U) 
+ 6n«'([/, Ux) + buU^U, U) + hsUxiU, U)^ 
+ hMU, U){U, Ux) + h.u^U, U)^ + heu{U, f/)3, 

■ hnUxxx + bisUUxxU + hgulU + b20UUxUx + h2iu'^Uxx 
+ h22U^UxU + hizu'^Ux + h^^u^U + h2^{U, U)Uxx + &26(f/, f/x)f/x 

+ b27{Ux, Ux)U + 628(t/, f/x.)f/ + h2,UUx{U, U)U 
+ 630«2^t/, U)Ux + &3lM2(f/, ^^^f; + h^2u\U, U)U 

+ 633(^7, [/)2f/, + 634(f/, f/)([/, + h^u^U, U)'U 

+ h,{U,U)^U, 

for which the following constraints guarantee the order to be 3 and the system not to be tri- 
angular: (61, 617) 7^ (0, 0) and at least one of feg, • • • , &16 and one of bis, • • • ? ^24, ^29, • • • , ^325 ^35 
must not vanish. However, when we consider a 3^*^ order symmetry for a 2^^^^ order system, 
we relax these constraints as follows (cf. section 12)): (&i,&i7) 7^ (0,0) and at least one of 



61, ... , biG and one of b 



17, 



636 must not vanish. 



Proposition 5.1. No 2"'' order system of the form j|l with a 3'"'^ order symmetry of the 
form or a 4*^* order symmetry exists. 
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Theorem 5.2. Any 3^'^ order system of the form (j^.jjj) with a S*'^ order symmetry has to 
coincide with either of the following two systems up to a scaling ofts, x, u, U {we omit the 
subscript oft^): 



' ut = {a + l){u^xx + 'iu^u^x + 'duul + Sw^Ma; + 3M^^(f/, U) 
+ 6«,(f/, Ux) + 3n,.(t/, Uf) + 2an(f/, U.,.,) 
+ (2a + 3)m(?7^, [/^.) + (10a + &)uxU^{U, U) + 2a'u3([/, U^) 
+ 6aM(f/, [/)([/, ?7^.) + au^U, U) + 2au3(f/, + a'u(f/, ?7)3, 

t/^ = f/... + 3([/, f/)t/,, + 6{U, U,)U, + 3{U,, U,)U + 3{U, UfU, 
— 2aUxxuU + (a + 3)v?JJ + GuUxUx + 3u'^Uxx — GaUxU^U 
+ 3u^Ux - 2aUxu{U, U)U - Aau^{U, Ux)U + 6u^{U, U)Ux 
-au^U - 2au^{U, U)U - au^{U, UfU, a: arbitrary, 

' ut = Uxxx + ^u^Uxx + ^uul + 3u'^Ux + 3Uxx{U, U) + Qux{U, Ux) 
+ 2u{U, Uxx) + MU.. Ux) + 10m,u2(?7, f/) + 2u3([/, [/^) 
+ 3ux{U, Uy + 6M(f/, U){U, Ux) + u^U, U) + 2u^{U, Uf 
' + u{U, U)\ 

Ut = -2uxxuU + ulU - GuxU^U - 2uxu{U, U)U - Au^{U, Ux)U 
-u^U -2u^{U,U)U -u'^{U,UyU 



Both system ()5.3|) and system ()5.4p admit the reduction U = 0. From this viewpoint, 
they are considered as generahzations of the Ibragimov-Shabat equation [36]. In addition, 
system fl5.3|) admits the reduction u = which changes it to a vector analogue of the 
Ibragimov-Shabat equation [15,34], 



Ut = Uxxx + 3(f/, U)Uxx + 6{U, Ux)Ux + 3{Ux, Ux)U + 3{U, UfUx. (5.5) 



We can linearize ()5.3p and ()5.4|) through the same change of dependent variables. In fact, 
both of them are 3'''^ order symmetries of a nontrivial 1*^* order system. 



which is naturally linearizable in the same way. 

5.2 Integrability of systems ( 15.31 ) and (15.41 ) 
5.2.1 System 

We note that system ()5.3|) possesses the following conservation law: 



{u' + {U,U)\ 

= [(a + l){2uuxx -ul + Qu'ux + m^) + 2{U, Uxx) - {Ux, Ux) 
+ (4a + Q)u^{U, Ux) + (2a + 6)m^m(C/, U) + (2a + 3)u^{U, U) 




(5.6) 



+ (a + 3)u\U, U)^ + 6(?7, U){U, Ux) + {U, [7)^] 



(5.7) 
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Then, if we introduce new variables w and W by 

yj = ^g/"(«2 + (C/,C/»dx' 



(5 i 



they satisfy a pair of hnear equations, 

Wt = {a + l)w. 



XXXI 

If we set = or M = 0, ()5.8p is reduced to the hnearizing transformation for the 
Ibragimov-Shabat equation [3,87] and that for its vector analogue ()5.5|) . respectively. 

5.2.2 System dOD 

System ()5.4|1 is obtained from ()5.3j) by rescaling t appropriately and taking the limit 
a ^ oo. As this fact implies in combination with ()5.7j) . system ()5.4|1 possesses the 
following conservation law: 

{u^ + {U, U))^ = {2uu^^ -ul + 6u\^ + u^ + 2u^u{U, U) 
+ Au^U, U,) + 2u\U, U) + Uf)^. 

Then, by the same change of variables as in section 15.2.11 

system (15. 4j) is decoupled into one linear equation and one trivial equation, 

= Wxxxi 

Wt = Q. 
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6 The case Ai = |, A2 = | — negative results — 



In this section, we search for 2 order and 3 order systems with a specific order sym- 
metry in the case of Ai = |, A2 = |. However, the resuhs turn out to be negative, as is 
shown below. 

The general ansatz for aAi = |, A2 = | homogeneous evolutionary system of 2°^^ order 
for a scalar function u and a vector function U takes the form 

r = aiu^^ + a2U^u^ + a^u'' + a4{U, U^) + a^^u^iU, U), 
{ Ut2 = aeUxx + a-jv^UxU + a^u^Ux + agu^U + aiou^{U, U)U. 

The following constraints guarantee the order to be 2 and the system not to be triangular: 

(ai, og) 7^ (0, 0), (a4, as) ^ (0, 0), (ay, og, Og, aio) ^ (0, 0, 0, 0). 

Similarly, the general ansatz for a 3'''^ order system takes the form 

' ut^ = hiUxxx + h2V?Uxx + hv?ul + hiU^Ux + hu^^ + hc{U, Uxx) 
+ h{Ux, Ux) + hu''ux{U, U) + hgu^U, Ux) + hou^{U, U) 
+ buu^{U,U)\ 

^ Ut3 = buUxxx + buU^UxxU + hi^uulU + hi^u^UxUx + bieU^Uxx (6-2) 
+ hnu^UxU + hsU^Ux + biQU^U + b2oUUx{U, U)U 
+ hiu\U, U)Ux + h2u\U, Ux)U + h2^u'>{U, U)U 
+ b2MU,U)^U, 

for which the following constraints guarantee the order to be 3 and the system not to be 
triangular: (61, 612) 7^ (0, 0) and at least one of b^, . . . , bu and one of 613, . . . , &24 must not 
vanish. However, when we consider a 3'''^ order symmetry for a 2"^^^ order system, we relax 
these constraints as follows (cf. section llj): (61, 612) 7^ (0, 0) and at least one of 61, ... , bu 
and one of 612, ... , &24 must not vanish. 

Proposition 6.1. No 2"^^ order system of the form ()6'. i|) with a 3'"'^ order symmetry of the 
form ()6'.i^ or a 4*^* order symmetry exists. 

Proposition 6.2. No 3''"' order system of the form ()6'.ijj) with a 5*'* order symmetry exists. 
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7 The case Ai = |, A2 = | 



In this section, we classify 2"^^^ order and 3'''^ order systems in the case of Ai = |, A2 = |. 
In the first part f section m|) . we present complete hsts of such systems with a specific 
order symmetry. In the second part (section 17. 2|) , we prove that the listed systems are 
linearizable. 



7.1 Lists of systems with a higher symmetry 

The general ansatz for a Ai = |, A2 = | homogeneous evolutionary system of 2'^'^ order 
for a scalar function u and a vector function U takes the form 

ut^ = aiu^^, + a2U^ + a^{U, U^^) + a4^{U^, f4) + a^u^^U, U) 

+ a^u^U, Uf + ajuiJJ. Uf + a^{U, U)\ (7.1) 
Up^ = flgf/^^ + aiou^U + auU^{U, U)U + ai2u{U, 11)^11 + ai3(f/, U)^U. 

The following constraints guarantee the order to be 2 and the system not to be triangular: 

(ai, 03, ag) ^ (0, 0, 0), (a3, 04, 05, ag, a?, ^s) (0, 0, 0, 0, 0, 0), 
(aio,aii,ai2) 7^ (0,0,0). 

Similarly, the general ansatz for a S''*^ order system takes the form 

' Ui3 = feiu^rxx + b2U^u.^ + h{U, U^^^) + hi{U^, U^^) + b^u'^u^iU, U) 

+ 6gM,(t/, Uf + 6iom(?7, f/)2(f/, [/,) + 6n(t/, f/)'(f/, t/x), 
^ Ut, = buU,,, + bnU^u.U + b^^u^U^ + bi^uu.,{U, U)U (7.2) 
+ 6i6M^(t/, f/)f/^ + bnu^{U, U^)U + bi8U^{U, UfU 

I +522(t/,f/)2(f/, 

for which the following constraints guarantee the order to be 3 and the system not to be 
triangular: (61, 63, 612) 7^ (0, 0, 0) and at least one of 63, ... , bn and one of 613, ... , 620 must 
not vanish. However, when we consider a 3^*^ order symmetry for a 2"^^ order system, we 
relax these constraints as follows (cf. section 12)): (&i, &3, &12) 7^ (0, 0, 0) and at least one of 
61, . . . , 611 and one of 612, ... , 622 must not vanish. 



Theorem 7.1. Any 2"*^ order system of the form ( |y.i| ) with a 3'"'^ order symmetry of the 
form ( |7.ijD has to coincide with the following system up to a scaling oft2, x, u,U {we omit 
the subscript 0/^2).' 

ut = u^. + 2(f/, U,^) + 2{U^, U^) + 2u{U, Uf + 2(f/, U)\ 
Ut = -u{U,UfU - {U,UfU. ^ ' ' 
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Theorem 7.2. Any 2"'^ order system of the form ( | 7. j| ) mt/i a 4*^^ order symmetry has to 
coincide with either l \7.!il^ or the following system up to a scaling oft2,x,u, U : 

^ut = -2{U, f/,.,.) - 2u^U, U) - 6u'{U, Uf - Qu{U, Uf - 2{U, U)\ 
\Ut = U^x + u^U + 3M2(f/, U)U + 3M(f/, U)^U + (f/, UfU. 



Theorem 7.3. Any S*"^ order system of the form ( |7.i4 ) 'W^^^ « S*^' or(/er symmetry has to 
coincide with either of the following two systems up to a scaling oft3,x,u, U {we omit the 
subscript oft^): 

Ut = u,,, + 2{U, U,,,) + 6{U,, U,,) + 2u,{U, Uf + 4(f/, Uf{U, U,), 
Ut = -u^{U, U)^U - 2{U, Uf{U, U,)U, 

Ut = u,,, + 2(f/, U^^^) + 6(f/,, f/,.,.) + 2n,.(f/, Uf + 4(f/, f/)3(t/, f/,), 

f/t = -n,.(f/, t/)2f/ - Au{U, U) {U, U^)U + Au{U, UfU, (7-6) 
-6([/,[/)2([/, + 4([/,[/)3[/,. 

We note that ()7.5|1 is the 3'''^ order symmetry of the 2"^^ order system ()7.3j) . 

7.2 Integrability of systems jH^^-iTrM 
7.2.1 Systems dUlD and (f73|> 

We present a procedure for solving system (j7.3p only, because its 3'''^ order symmetry (j7.5j) 
can be solved in the same way. For system (j7.3j) . if we introduce a new variable w by 

w = u+{U,U), (7.7) 

it solves the linear equation, 

Wt = Wxx- 

Once we know w{x, t) by solving this equation, we obtain from the relation (([/, f/)~^)^ = 
Aw that 

1 /■* 1 

4 / w{x,t')dt' + 



{U{x,t),U{x,t)f 7o (f/(a;,0),f/(a;,0))2- 

Then we can determine u{x^ t) by using ()7.7|) . Finally, noting the relation ((t/, U)^^U)^ 
0, we obtain the following expression for U{x,t): 

U{x,t) = tU{x,0). 

1 + 4:{U{x, 0), U{x, 0))2 w{x, t')dt'] ' 
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7.2.2 System (fTTH) 

For system ()7.4|) . we have the relation {u + {U, U))t = 0. Thus, we can set 

u+{U,U)=(P{x), 

where the function does not depend on t. Then the equation for U is rewritten in 
terms of 

Ut = f/.. + (p^U. (7.8) 

The solutions of ()7.8|) are given by 

U{x,t) = j dAe^*^(a;; A), 

where A) is a solution of the ordinary differential equation 

+ 0=^^ = A^. 

The following commutation relation indicates that system ()7.4j) possesses a polynomial 
higher symmetry of every even order (cf. [7,10]): 

[dl + ^\ (9,2 + 0^)"] = 0, n = l,2,... . 

7.2.3 System ^TM 

For system ()7.6|) . if we introduce a new variable w by 

w = u+ {U,U), 

it solves the linear equation. 

Once we know w{x,t), we obtain from the relation {{U, f/)^^)^ = Awx that 

{U{x,t),U{x,t))^ Vo "'^■^■^'^^^^ ■ ([/(x,0),[/(x,0))2 
Then the equation for {U,U)^^U can be rewritten as 

^ =^u+{U,U)){U,U)'(^=L=U 



4 / ti;,(x,t')dt' + 777— -^.r— (7.9) 



4w;(a;,t) / 1 



U 



The general solution of this equation is given by 

^ -.U, = f, (a [ w{x,t')dt' + [ .^--^—-^dx' 



^WU) V io J (f/(x',0),f/(a;',0))2 

j = l,2,...,iV, (7.10) 
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where fi{z), . . . ^ f]\i{z) are arbitrary functions of z, except that they must satisfy one 
constraint, Ej=i [/j(-s)] = 1- Combining (fTTTH) with ((7H), we arrive at the following 
formula: 

Ujix,t) = -^f,iO. J = l,2,...,iV, 
where ^{x, t) = A w{x, t')dt' + r{U{x', 0), U{x', 0))~W. 
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8 Concluding remarks 



In this paper, we have presented a classification of integrable evolutionary systems in 1 + 1 
dimensions for one scalar unknown u{x,t) and one vector unknown U{x,t). We focused 
on polynomial systems that are homogeneous under a suitable weighting of dx, dt, u{x, t), 
U{x,t), and considered five distinct weightings for u, U relative to a fixed weight of d^- 
Then, with the help of a computer algebra program, we obtained the complete lists, up to 
a scaling of variables, of 2°^^ order systems with a 3'"'^ order or a 4**^ order symmetry and 
3rd Qj-(^Qj- systems with a 5^^ order symmetry. We demonstrated the integrability of nearly 
all listed systems by constructing a Lax representation or a linearizing transformation, 
or, in some cases, by identifying an integrable closed subsystem contained in the system 
under investigation. 

The following tables give a quick overview of the systems found. Note we use "MT" as 
an abbreviation for "Miura-type transformation" , including Miura map plus potentiation. 
In these tables, we set the weight of dx at unity, without any loss of generality. For 
full details regarding Lax representations, transformations, references, etc., the reader is 
referred to the corresponding part of the paper identified through the equation number. 
Here, we would like to make a few remarks on our classification results: 

• The most interesting classification results are obtained for the case Ai = A2 = 1, 
namely, the Burgers/pKdV/mKdV weighting. The lists in this case consist of a 
large number of systems, which are shown to have a very wide variety of underlying 
structures. We compared these lists thoroughly with the lists of two-component sys- 
tems by Foursov-Olver [18,25,26], refined and generalized their work, as described 
in the introduction. 

• We found a number of pairs/triplets of scalar- vector systems connected through 
transformations of dependent variables. Besides standard Miura transformations 
that map both the scalar and vector variables to new ones (see, e.g. ()4.44j) ). we 
also found Miura-type transformations that act only on the scalar variable and do 
not change the vector variable (see, e.g. ()4.46|) combined with potentiation v = Ux). 
For some other systems, we showed that a new scalar variable defined in terms of 
the old scalar and vector variables satisfies a closed integrable equation, such as the 
KdV equation or a linear equation. 

• The search for such transformations in our case of scalar-vector systems is simple 
in comparison to scalar-scalar systems. For instance, the ansatz that a new scalar 
variable depends on the original vector variable U only through scalar products 
{dl^U, d'^U) narrows down the candidates for such transformations considerably. 
This leads us to the counter-intuitive observation that scalar-vector systems are, 
in a sense, more tractable than scalar-scalar systems. This is probably one reason 
why, unlike our work, Foursov and Olver proved integrability^^ for only a small 
proportion of their two-component systems [18,25,26]. 

Finally, we mention some problems not solved in this paper: 
^^They discussed the existence of a recursion operator or a bi-Hamiltonian structure. 
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• How can the integrability of the three systems (j4.5j) . (j4.8p and (j4.9|) be estabhshed 
along the hnes of this paper? The main obstacle is that we know neither a linearizing 
transformation nor a proper Lax representation for the two-component Burgers 
system ()4.35|) . The dependence of the functional form of travelling- wave solutions 
on the boundary conditions and the velocity implies that (j4.35|) is a highly nontrivial 
system and not linearizable by a naive extension of the Hopf-Cole transformation. 

• Some scalar-vector systems are converted to a triangular form, i.e. a closed subsys- 
tem plus remaining equations coupled to it. When the remaining equations contain 
a nonlinear PDE in its own variable (cf. ()4.54b|) or ()4.58b|) ). it seems to be especially 
difficult to solve them explicitly for a given solution of the subsystem. Is there any 
method, like an extension of the inverse scattering method, for dealing with such 
triangular systems analytically? 

• Can one construct an explicit formula for the general solution of ()4.55|) ? This 
equation is obtained from system ()4.20j) or 1)4.211) by converting it to a triangular 
form and then considering the special case in which the solution of the subsystem, 
KdV equation in this case, is identically zero. 

Although we concentrated our attention on the five distinct weightings for u, U in 
this paper, we also found integrable systems that are homogeneous under a different 
weighting of variables. Namely, we obtained systems of coupled KdV-mKdV type, e.g. 
fl07|l . (lOnjl and (PT7|1 ^^ together with the proof of their integrability. We are 

planning to complete a classification of integrable systems of this type, i.e. scalar-vector 
systems with weights Ai = 2, A2 = 1, in a subsequent paper. Preliminary results can be 
viewed on the web page 

http : //lie . math . brocku . ca/ twolf /htdocs/sv/over . html . 



^^We note that, in addition to a scaling of variables, these systems admit another equivalence trans- 
formation u = u + k{U, U), U = U. 
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Table 4: An overview of the considered classes with unit weighting of d„ 



weights 
(Ai.Aa) 
of u, U 


weights of 
dt, dr in 
sys., sym. 


system 


comments 


(2,2) 


2, 3 


none 




2, 4 


none 




3, 5 


< 


-«t = {U,U^), 

, C/t = + UxU + 2uUx- 


multi-component generalization of a 
Drinfel'd— Sokolov system [43,44], see 
[45-47] 


mt < 


' Ut = Uxxx + 6uUx — 6(U, Ux), 

, Ut = Uxxx + duxU + euUx- 


known sls 3, Jord&n KdV system [27,28, 
34,48] 


< 


' Ut = Uxxx + 3uUx + 3(11, Ux), 

, Ut = UxU + uUx- 


multi-component generalization of 
Zakharov— Ito system [52,53], see [54] 




' Ut = Uxxx + 6uUx — 12{U, Ux), 

. Ut = —2Uxxx — GuUx- 


multi-component generalization of 
Hirota-Satsuma system [57], see [58] 


(1,1) 


2, 3 
or 
2, 4 


(m < 


' Ut = ^{l + 2a){uxx + 2uUx) + ^{U,Ux), 

Ut = Uxx + l{l- a)uxU + uUx + 1^(1 - 4a)u'^U 
- i(C/, U)U, a : arbitrary. 


linearized by an extended Hopf— Cole 
transformation 




' Ut = Uxx + 2uux + 2{U, Ux), 

^ Ut = -\uxU - \u'^U - \(U, U)U. 


is scaling limit of 14. 3i. linearized by 
the same transformation 


m < 


' Ut = Uxx + 2uUx + ((7, Ux), 

Ut = \uxU + uUx- 


contains two-component Burgers sys- 
tem 14.351 as closed subsystem, inte- 
grability unproven 


3, 5 




' Ut = a{uxxx + Suuxx + 3ul + Su'^Ux) + Ux{U, U) + 2u(U, Ux) 
+ 2{U,Uxx) + 2{Ux,Ux), 
Ut = Uxxx + |(1 - a)uxxU+ luxUx + \uUxx + |(1 - 2a)uuxU 
+ lu^Ux - (U, Ux)U + 1(1 - Aa)u^U - ^u{U, U)U, 

a : arbitrary. 


is symmetry of 14. 3i 


(HI < 


■ Ut = Uxxx + 3uUxx + 3Ux + 3u^Ux + Ux{U, U) + 2u{U, Ux) 

+ 2(U,Uxx) + 2(Ux,Ux), 
. Ut = -\uxxU - luuxU - (U, Ux)U - \u-^U - \u{U, U)U. 


is symmetry of 14. 4i. scaling limit of 


.—1 J 

lot < 


' Ut = Uxxx + 3uUxx + 3u1 + 'iu^Ux + Ux{U, U) + 2u(U, Ux) 

1 / T T T T \ 1 It T T T \ 

+ {U, Uxx) + [Ux, Ux), 

. Ut = \uxxU + UxUx + uuxU + u'^Ux + ^{U, U)Ux + ^{U, Ux)U. 


is symmetry of 14. 5i. see there 


<oi < 


' Ut = Uxxx + SuUxx + 3u'^ + 3u^Ux + Ux{U, U) + 2u{U, Ux) 

+ {U,Uxx) + {Ux,Ux), 
. Ut = \uxxU + UxUx + UUxU + U^Ux + {U, U)Ux- 


contains 3"^ order symmetry of two- 
component Burgers system I4.35i. in- 
tegrability unproven 




' Ut = 3ux{U,U) + 3{U,Uxx) - 3{U,U)'^ , 

. Ut = Uxxx + UxxU + UxUx- 3{U, Ux)U. 


obtamed irom 14.121 by iVl 1 




' Ut = 2ux{U, U) + 2(U, Uxx) - {Ux, Ux) - 2{U, U)^ , 
. Ut = Uxxx + UxxU + 2uxUx - 2{U, U)Ux - 2{U, Ux)U. 


linearizable by change of variable, re- 
lated to Kaup-Kupershmidt equation 
in a certain manner 


i4J2t 


' Ut = Ux(U, U) + 2u(U, Ux) + {U, Uxx) + {Ux, Ux), 

. Ut = Uxxx + UxxU + UxUx- 2uUxU - U^Ux + (U, U)Ux - (U, Ux)U. 


connected witn lo.ol ana I4.1UI by M i 




' Ut = Uxxx + \ul + ^{Ux,Ux), 
Ut = UxUx- 


is potential form of 13.51 




' Ut = Uxxx + 3ul + 2aux{U, U) + a{U, Uxx) + a(Ux,Ux) + b{U, U)^, 
. Ut = UxxU + 2uxUx + a{U, U)Ux + a{U, Ux)U, (a, b) ^ (0, 0). 


for 6 7^ a^/4, transformed to 13. 5i: for 
b = a^/4, to KdV equation + linear 
vector equation coupled to it 
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Table 5: Continuation 



weights 
of u, U 


weights of 
dt , O'r in 
sys., sym. 


system 


comments 


(1, 1) 


3, 5 


KWl . 


' ut = Uxxx + 3n^ - 'i{Ux, Ux), 

, Ci = Uxxx + GUxUx- 


is potential form of 13. 4i. MT connects 
to iTITl 








' Ut = Uxxx + + Ux{U, U) + (U, Uxx), 
. Ut = Uxxx + SUxxU + SUxUx + {U, Ux)U. 


obtained from lOSl by MT 






Bnn • 


' Ut = Uxxx + 3ul + 2Ux{U, U) + {U, Uxx) + ^{Ux,Ux), 

Ut = Uxxx + GuxxU + GuxUx + 2(1/, Ux)U. 


obtained from by MT 






■ 


Ut = Uxxx + 3ul + AUx(U, U) + 2{U, Uxx) + {Ux, Ux) + 

Ut = -2Uxxx - GuxxU - GuxUx - MU, Ux)U. 




obtained from lOUl by MT 






BHTH ■ 


' Ut = Uxxx + ul- 12{U, Uxx) + l2(Ux,Ux) - 4(C/, (7)2, 
^ Ut = iUxxx + UxxU + 2uxUx + 4(1/, U)Ux + 4(1/, Ux)U. 


reduced to triangular system: KdV 
equation + linear vector equation cou- 
pled to it 








' Ut = Uxxx - \u^Ux + ^Ux{U,U) +u{U,Ux) + {U,Uxx) 

Ut = -UxUx - ^u^Ux + l{U,U)Ux. 


+ {Ux,Ux), 


converted to triangular system: KdV 
eq. + nonlinear eq. with interesting re- 
duction 14.551 -I- linear vector eq. 






ion . 


' Ut = Uxxx - \u^Ux + lux{U,U) +u{U,Ux) + (U,Uxx) 

. Ut = -UxUx - \u^Ux + \{U,U)Ux + {U,Ux)U. 


+ {Ux,Ux), 


exactly like for i4.20i only different lin- 
ear vector eq. 






1 


' Ut = Uxxx - lu'^ux + \ux{U,U) +u{U,Ux) + {U,Uxx) 
Ut = UxxU + UxUx - UUxU - \u^Ux + \ {U,U)Ux 

+ {U,Ux)U. 


+ {Ux,Ux), 


obtained in [69], converted to trian- 
gular system: KdV equation + linear 
equations coupled to it, admits defor- 
mation connected with 13. 5i by MT 






to^ ■ 


' Ut = Uxxx - lu^Ux + lux{U, U) + U{U, Ux) + (C/, Uxx) 

+ ^{U,U)\ 

. Ut = -UxUx - ^u^Ux - i(C/, U)Ux + ^u{U, U)U. 


+ {Ux,Ux) 


contains interesting triangular system 
I4.58i as closed subsystem: KdV equa- 
tion + nonlinear equation coupled to it 






ion . 


' Ut = Uxxx - lu^Ux +Ux{U,U) +u{U,Ux) + {U,Uxx) 4 

- luHu,U) + ^{U,U)\ 
. Ut = ^UxxU + ^{U,Ux)U - ^u^U + ^u{U,U)U. 


-(Ux,Ux) 


is symmetry of l^'^ order system I4.61i. 
converted to triangular system: KdV 
eq. -j- Riccati eq. coupled to it -\- linear 
vector eq. coupled to them 






iTIEl < 


' Ut = Uxxx + U^Ux + Ux{U, U), 
.Ut = Uxxx+U^Ux + {U,U)Ux. 


equivalent to a single vector mKdV 
equation for vector [u, U) 






KM < 


' Ut = Uxxx + 2u^Ux + Ux(U, U) + u{U, Ux), 

. Ut = Uxxx + uuxU + u^Ux + {U, U)Ux + {U, Ux)U. 


equivalent to a single vector mKdV 
equation for vector (u, U) 






iOrl . 


' Ut = Uxxx — Gu^Ux + 6Ux{U, U) + l2u{U, Ux), 

, Ut = Uxxx - 12uuxU - Gu^Ux + 6{U, U)Ux. 


known as a Jordan mKdV system [27] , 
connected with and BTT^ by MT 








' Ut = Uxxx - Gu'^Ux + Ux{U, U) + 2u{U, Ux) + {U, Uxx) + {Ux,Ux), 

Ut = Uxxx + 3uxxU + SuxUx - duuxU - 3u'^Ux + {U, U)Ux 
+ 3{U,Ux)U. 


admits Lax representation, connected 
with BTTSl by MT 






iTM ■ 


' Ut = Uxxx - Gu^Ux + Ux{U, U) + 2u{U, Ux) + (U, Uxx) + (Ux,Ux), 

Ut = Uxxx + GuxxU + GuxUx - l2uuxU - Gu^Ux + {U, U)Ux 
+ 'l{U,Ux)U. 


multi-component generalization of a 
modified Jaulent— Miodek flow [83] , ad- 
mits Lax representation, connected 
with Kvfi by MT 






ton ■ 


' Ut = Uxxx - &U^Ux +Ux{U,U) + 2u{U,Ux) + {U,Uxx) - 

Ut = —2Uxxx — duxxU - GuxUx + 12uuxU + Gu^Ux + 
-2{U,Ux)U. 


¥(Ux,Ux), 
(U,U)Ux 


connected with and BTT^ by MT 



53 



Table 6: Continuation 



weights 
(Ai,A2) 
of u, U 


weights of 
dt, Or in 
sys., sym. 


system 


comments 




2, 3 


none 






2, 4 


none 






3 5 


■ 




= (a + \){Uxxx + "iu^Uxx + 9UM^ + "iu^Ux + 3na;2!(r/, 

+ %u^{U, U^) + 3w^(r/, C/>2) + 2au{U, Ua:^) 

+ (2a + ?,)u{Ux, Ux) + (10a + 6)«^n2((7, !7) + 2a«3(;7, f/^) 

+ 6aM(C/, 6f)(C/, Lfo;) + au^{U, U) + 2au^U, UY + a.u{U, U)^ , 

= Ua:xx + 3(C/, U)Ua:^ + &{U , Ua:)Ua: + 3{Ua:, U^)U + 3{U, U)^U^ 

— 2aUxxuU + (a + 3)u^U + GuUxUx + Su^Uxx — GaUxU'^U 
+ 3u^Ux - 2auxu{U, U)U - Aau'^{U, Ux)U + Gu'^{U, U)Ux 

- au^U - 2au^{U, U)U - au'^{U, U)^U, a : arbitrary. 


is symmetry of 1^^ order system ^^^J , 
extension of vector Ibragimov— Shabat 
equation, linearizable by change of 
variables 






El ■ 


Ut 


= Uxxx + Su'^Uxx + Quul + 3u^Ux + 3Uxx{U, U) + 6Ux{U, Ux) 

+ 2u(U,Uxx) + 2u{Ux,Ux) + 10uxu'^{U,U) + 2u^(U,Ux) 
+ 3ux(U,U)^ + 6u{U,U){U,Ux) + u^(U,U) + 2u^{U,U)^ 
+ u{U,U)^, 

= -2uxxuU + ulU - GuxU^U - 2uxu{U, U)U - 4u^{U, Ux)U 
-u^U - 2u^{U, U)U - u^(U, U)^U. 


is symmetry of f*' order system 15.61. 
scaling hmit of 15. 3i. linearized by the 
same change of variables 




2, 3 


none 






2, 4 


none 






3, 5 


none 






2, 3 


O • 


. Ut 


= uxx + 2(U, Uxx) + 2(Ux, Ux) + 2u{U, C/>3 + 2{U, U)^ , 
= -u(U,U)^U- (U,U)''^U. 


ultralocal change of variables gives lin- 
ear equations 












2, 4 


El • 


. Ut 


= -2{U, Uxx) - 2ui{U, U) - 6u'2{U, (7)2 - 6u{U, U)^ - 2{U, C/)*, 
= Uxx + u^U + 3u^{U, U)U + 3u(U, UfU + (U, UfU. 


ultralocal change of variables gives lin- 
ear equations 




3, 5 


ESJ • 


' Ut 

. Ut 


= Uxxx + 2(U, Uxxx) + e(Ux, Uxx) + 2Ux(U, C/>3 + 4(U, U)^U, Ux), 

= -Ux{U,U)^U -2{U,U)^{U,Ux)U. 


is symmetry of 17.31 






EU ■ 


' Ut 

I'" 


= Uxxx + 2(U, Uxxx) + 6{Ux, Uxx) + 2ux{U, C/>3 + 4(U, U)^{U, Ux), 

= -Ux(U, U)^U - 4u(!7, U){U, Ux)U + 4u{U, U)^Ux 
- 6(U, U)^{U, Ux)U + A{U, U)^Ux. 


ultralocal change of variables gives lin- 
ear equations 
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